DEPARTMENT  OF  THE  AIR  FORCE 
AIR  UNIVERSITY 

AIR  FORCE  INSTITUTE  OF  TECHNOLOGY 


Wright-Patterson  Air  Force  Base,  Ohi 


4 


4 


AFIT/GE/ENG/90D-19 


NUMERICAL  APPPROXIMATIONS 
OF  THE  GREEN’S  FUNCTIONS 
FOR  A  THEORETICAL  MODEL 
OF  AN  APERTURE  FED  STACKED-PATCH 
MICROSTRIP  ANTENNA 

THESIS 

Ronald  E.  Erwert 
Captain,  USAF 

AFIT/GE/ENG/90D-19 


Approved  for  public  release;  distribution  unlimited 


AFIT/GE/ENG/90D-19 


NUMERICAL  APPROXIMATIONS  OF  THE  GREEN’S  FUNCTIONS 
FOR  A  THEORETICAL  MODEL  OF  AN  APERTURE  FED 
STACKED-PATCH  MICROSTRIP  ANTENNA 


THESIS 


Presented  to  the  Faculty  of  the  School  of  Engineering 
of  the  Air  Force  Institute  of  Technology 
Air  University 
In  Partial  Fulfillment  of  the 
Requirements  for  the  Degree  of 
Master  of  Science  in  Electrical  Engineering 


Ronald  E.  Erwert,  B.S. 
Captain,  USAF 


December,  1990 


Approved  for  public  release;  distribution  unlimited 


Preface 


The  purpose  of  this  study  was  to  develop  numerical  approximations  for  the  Green’s  functions 
of  a  theoretical  model  of  an  aperture  fed  stacked-patch  microstrip  antenna.  Each  complex  valued 
Green’s  functions  was  approximated  by  three  functional  forms.  One  component  of  each  Green’s 
function  was  well  approximated  by  a  single  function.  The  other  component  required  a  piece- 
wise  approximation  but  with  each  approximating  function  valid  over  the  same  range,  of  radial 
separation,  for  each  function.  All  the  approximations  were  determined  using  a  least  squares  curve 
fitting  method  and  excellent  agreement  with  the  original  data  was  obtained.  These  approximations 
were  given  to  Capt  Leon  Irvin  for  use  in  a  separate  effort  to  determine  the  patch  surface  currents 
using  a  method  of  moments  analysis.  Approximations  of  this  type  should  provide  for  a  much  more 
computationally  efficient  method  of  moments  analysis. 

I  would  like  to  thank  my  advisor.  Major  Harry  Barksdale,  for  his  help,  encouragement,  and 
patience  throughout  the  last  year.  I  would  also  like  to  thank  the  members  of  my  thesis  committee, 
Capt  Philip  Joseph  and  Capt  Byron  Welsh  for  their  help  as  well  as  patience  and  understanding. 

The  last  year  and  a  half  have  been  some  of  the  best  and  worst  times  of  my  life.  The  association 
with  the  members  of  the  LO  track  will  stand  as  a  highlight  of  my  life.  In  what  could  have  been  a 
competitive  race  for  position  in  our  classes,  the  members  of  the  LO  track  were  always  ready  and 
willing  to  help  and  encourage  each  other.  The  entire  Air  Force  could  benefit  from  this  cooperative 
type  of  relationship.  On  the  downside,  the  last  year  and  a  half  have  been  the  worst  in  regards 
to  personal  and  family  life.  I  would  like  to  give  special  thanks  to  my  wife  Carole,  for  her  love, 
patience,  and  understanding  over  the  last  year  and  a  half.  She  sacrificed  her  career  that  I  might 
obtain  this  degree  and  words  will  never  be  able  to  express  iny  appreciation.  I  Would  also  like  to 
give  thanks  to  my  sons  David  and  Jonathan  for  doing  without  “Daddy”  much  of  the  time. 

In  any  undertaking,  such  as  a  thesis,  hind  site  is  always  much  clearer  than  foresight.  At 
the  outset  of  this  effort  an  unwritten  goal  was  to  rewrite  the  code  used  by  Nazar  to  eliminate  the 


n 


use  of  commercial  software  packages.  The  primary  goal  was  determining  the  Green’s  functions 
approximations.  Over  the  summer  a  great  deal  of  time  was  spent  investigating  different  analysis 
and  integration  techniques.  During  this  time  I  seemed  to  lose  sight  of  the  primary  goal,  obtaining 
the  approximations.  When  the  decision  was  made  to  press  on,  using  the  routines  developed  by 
Nazar,  much  valuable  time  had  been  lost.  A  clearer  plan  of  attack,  with  a  sequence  of  important 
events,  might  have  kept  me  on  track  better.  The  fault  rests  with  the  author  and  not  the  advisor. 
Sometimes  the  lessons  hardest  learned  are  the  best  remembered. 


Ronald  E.  Erwert 
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AFIT/GE/ENG/90D-19 

Abstract 

The  numerical  approximations  of  the  Green’s  functions  for  a  theoretical  model  of  an  aperture 
fed  stacked^tch  microstrip  antenna  are  presented.  The  Green’s  functions  approximated  are  those 
developed  by  Nazar  (AFIT  Thesis  GE-89D-37),  '‘Green’s  Functions  for  a  Theoretical  Model  of 
an  Aperture  Fed  Stacked-Patch  Microstrip  Antenna”.  Each  complex  valued  Green’s  function  was 
written  as  a  function  of  the  source-observer  separation.  The  behavior  of  each  Green’s  function  was 
analyzed,  and  it  was  found  that  they  could  all  be  approximated  by  three  separate  functions.  The 
approximating  functions  were  either  polynomials  or  rational  functions  composed  of  polynomials. 

One  component  could  be  approximated  by  a  single  function  over  the  entire  range  of  radial 
separations,  while  the  other  required  a  piecewise  approximation  over  the  same  intervals  for  each 
function  approximated.  The  coefficients  of  the  functions  were  determined  by  using  a  least  squares 
curve  fitting  method  and  were  used  in  a  separate  effort^AFIT-Theeis-G  E.9Q]3::2j^  to  perform 
a  moment  method  analysis  of  the  antenna.  Custom  written  FORTRAN  code  was  developed  to 
determine  the  coefficients  of  the  approximations.  j 

Example  results  for  several  of  the  Green’s  functions  approximated  are  presented^The  agree¬ 
ment  between  the  approximations  and  the  original  functions  was  excellent. 
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NUMERICAL  APPROXIMATIONS  OF  THE  GREEN'S  FUNCTIONS 
FOR  A  THEORETICAL  MODEL  OF  AN  APERTURE  FED 
STACKED-PATCH  MICROSTRIP  ANTENNA 


/.  Introduction 

The  microstrip  antenna  concept  was  first  proposed  by  Deschamps  in  1953  [1].  Tlie  concept 
lay  dormant  until  the  early  1970’s  when  there  was  an  immediate  demand  for  low-profile  antennas 
on  a  new  generation  of  missiles  [2]  and  better  theoretical  models  and  manufacturing  techniques 
were  developed.  Low  cost,  ease  of  construction,  thin  profile,  and  modular  design  are  among  tlie 
numerous  advantages  offered  by  microstrip  antennas.  Conformability  makes  them  easy  to  integrate 
into  the  skin  of  an  aircraft  or  missile  without  degrading  the  vehicle  aerodynamics.  E.xtensive  efforts 
continue  to  improve  the  designs,  integrate  them  into  new  applications,  and  develop  better  analysis 
techniques  to  predict  antenna  performance 

J.l  Microsirip  Anienna  Background 

1.1.1  Definition.  The  basic  geometry  of  a  microstrip  antenna  consists  of  a  conducting  strip 
radiator  or  patch  printed  on  a  grounded  dielectric  substrate  (see  Figure  1).  The  antenna  input  or 
feed  can  be  either  a  coaxial  cable  or  a  microstrip.  The  antenna  patch  conductors  can  in  principle  be 
any  shape  but  are  usually  rectangular  or  circular  to  simplify  analysis  and  performance  prediction. 
Generally  the  substrate  dielectric  constant  should  be  low  (e  «  2.5)  to  enhance  the  fringe  fields 
which  account  for  the  radiation  [1]. 


Antenna  Patch 


Antenna  Patch 


Figure  1.  Microstrip  Antennas  with  (a)  Coaxial  Feed  and  (b)  Microstrip  Feed  [3]. 

l.l.S  Advantages.  Compared  to  conventional  antennas,  microstrip  antennas  have  numerous 
advantages  and  can  be  used  in  a  wide  variety  of  applications  over  a  frequency  range  of  approximately 
100  MHz  to  50  GHz.  Some  of  the  major  advantages  are  [1]: 


•  light  weight,  low  volume,  low  profile  planar  configurations  which  can  be  made  conformal 

•  low  fabrication  cost,  readily  amenable  to  mass  production 

•  can  be  made  thin;  hence  they  do  not  perturb  the  aerodynamics  of  the  host  aerospace  vehicles 

•  the  antennas  may  be  easily  mounted  on  missiles,  rockets  and  satellites  without  major  alter¬ 
ations 

•  microstrip  antennas  are  compatible  with  modular  designs  (Solid  state  devices  such  as  oscil¬ 
lators,  amplifiers,  variable  attenuators,  switches,  modulators,  mixers,  phase  shifters,  etc.  can 
be  added  directly  to  the  antenna  substrate  board) 

•  feed  lines  and  matching  networks  may  be  fabricated  simultaneously  with  the  antenna 


1.1.3  Disadvantages.  Unfortunately  microstrip  antennas  do  have  some  disadvantages  com' 
pared  to  conventional  antennas  such  as  [1]: 


•  narrow  bandwidth 

•  loss,  hence  somewhat  lower  gain 

•  most  microstrip  antennas  radiate  into  a  half  plane 

•  poor  end-fire  radiation  performance 

•  poor  isolation  between  the  feed  and  radiating  elements 
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1.2  Problem  Background 


As  previously  stated,  narrow  bandwidth  is  one  of  the  disadvantages  of  microstrip  antennas. 
The  bandwidth  of  the  antennas  shown  in  Figure  1  is  typically  1-5%  of  the  resonant  frequency. 
Increasing  the  substrate  thickness  can  increase  the  bandwidth  but  any  bandwidth  gains  are  offset 
by  loss,  radiation,  and  impedance  mismatch  problems  [3]. 

Captain  James  Nazar  developed  an  analytical  model  of  a  configuration  combining  two  band¬ 
width  improving  techniques,  aperture  feeding  from  a  microstrip  line  and  using  stacked  patch  radi¬ 
ating  antenna  elements.  A  picture  of  this  aperture  fed  stacked-patch  microstrip  antenna  is  shown 
in  Figure  2.  In  his  thesis  he  presented  a  theoretical  model  for  this  antenna  using  the  mixed  po¬ 
tential  integral  equation  (MPIE)  approach.  The  Green’s  functions  associated  with  the  vector  and 
scalar  potentials  are  evaluated  in  the  spatial  domain  using  stratified  media  theory.  He  outlines  a 
method  of  moments  technique  to  solve  for  the  currents  but  no  actual  solutions  for  the  currents  are 
calculated. 

The  Green’s  functions  developed  are  expressed  in  the  form  of  Sommerfeld  integrals.  The 
moment  method  analysis  requires  surface  integration  of  the  Green’s  functions  over  the  source  and 
observer  coordinates.  The  Green’s  function  integrands  are  dependent  on  the  radial  separation 
between  the  source  and  observer.  To  numerically  evaluate  the  integrals  of  each  Green’s  function 
at  every  value  of  radial  separation  required  in  the  surface  integrations  would  be  extremely  time 
consuming.  Nazar  proposed  numerically  evaluating  the  integral  Green’s  functions  over  a  finite 
number  of  radial  separations  and  interpolating  to  find  the  intermediate  values.  He  states  polynomial 
averaging  may  be  an  acceptable  interpolation  technique  [3]. 

l.S  Problem  Siaiement 

In  this  thesis,  a  subset  of  the  Green’s  functions  developed  by  Nazar  will  be  numerically 
evaluated  at  a  finite  number  of  observer-source  radial  separations  and  approximations  for  each 
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1.5  Research  Questions 


Answers  to  the  following  questions  will  be  found: 

1.  What  are  the  form  functions  necessary  to  accurately  model  the  Green’s  functions? 

2.  How  well  do  these  functions  approximate  the  original  Green’s  function? 

1.6  Assumptions 

Because  this  thesis  effort  is  a  continuation  of  Nazar’s  work,  the  same  simplifying  assumptions 
were  made  about  the  material  parameters  to  make  the  mathematical  analysis  tractable.  Each 
dielectric  is  considered  isotropic,  homogeneous,  and  lossless.  The  dielectric  layers  and  ground  plane 
are  considered  infinite  sheets.  The  ground  plane  and  antenna  patch  conductors  are  infinitely  thin, 
perfect  electrical  conductors,  while  the  dielectric  layers  have  finite  thickn^s.  Finally,  the  Green’s 
functions  derived  by  Nazar  were  assumed  correct. 

1. 7  Scope  and  Limitations 

Only  the  aperture  fed  stacked-patch  microstrip  antenna  was  analyzed.  The  unique  Green’s 
functions  derived  by  Nazar  were  evaluated  numerically  and  approximated  by  appropriate  polyno¬ 
mials.  Once  these  functions  were  determined,  the  coefficients  were  given  to  Capt  Irvin  for  use  in  a 
method  of  moments  analysis  to  determine  the  electric  and  magnetic  currents  of  the  antenna. 

1.8  Thesis  Organization 

The  rest  of  this  document  is  organized  in  the  following  manner.  Chapter  II  is  a  review  of 
several  current  articles  on  aperture  fed  microstrip  antennas,  stacked-patch  microstrip  antennas, 
aperture  fed  stacked-patch  microstrip  antennas,  and  the  method  of  full  wave  analysis.  Chapter 
III  contains  the  theory  necessary  to  develop  the  Green’s  functions  in  the  theoretical  model  used 
to  analyze  this  antenna.  Chapter  IV  presents  an  explanation  of  what  functions  are  approximated. 
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what  the  approximations  are  and  how  well  the  functions  approximate  the  Green’s  functions.  The 
conclusions  and  recommendations  are  discussed  in  Chapter  V.  Finally  the  components  of  the  vector 
potential  derived  in  Chapter  III  as  well  a  complete  list  of  the  Green’s  functions  used  in  the  model 
are  contained  in  appendices. 
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II.  Literature  Review 


This  literature  review  summarizes  current  work  on  microstrip  antennas.  Since  this  thesis 
effort  concentrates  on  further  development  of  a  theoretical  model  for  an  aperture  fed  stacked-patch 
microstrip  antenna,  only  current  reports  on  aperture  fed  microstrip  antennas,  stacked  patch  mi¬ 
crostrip  antennas,  aperture  fed  stacked-patch  microstrip  antennas,  and  full  wave  analysis  techniques 
will  be  summarized.  The  final  work  reviewed  is  Captain  James  Nazar’s  thesis,  the  departure  point 
for  this  thesis  effort. 


S.J  Aperiure  Coupled  Antennas 

In  1985  D.  M.  Pozar  [6]  presented  a  new  technique  for  feeding  microstrip  antennas.  He 
proposed  coupling  a  microstrip  patch  on  a  substrate  to  a  microstrip  line  feed  located  on  a  parallel 
substrate  via  an  aperture  in  the  ground  plane  separating  the  two  substrates  (see  Figure  3).  To 
make  the  feedline  location  easier  to  visualize,  the  feed  substrate  in  Figure.3  is  shown  transparent. 
Pozar  cited  three  advantages  of  an  aperture  coupled  configuration: 


1.  The  configuration  is  well  suited  for  monolithic  phased  arrays,  where  active  devices  can  be 
integrated  on,  for  example,  a  gallium  arsenide  substrate  with  the  feed  network,  and  the  radi¬ 
ating  elements  can  be  located  on  an  adjacent  (low-dielectric  constant)  substrate,  and  coupled 
to  the  feed  network  through  apertures  in  the  ground  plane  separating  the  two  substrates.  The 
use  of  two  substrates  thus  avoids  the  deleterious  effect  of  a  high-dielectric-constant  substrate 
on  the  bandwidth  and  scan  performance  of  a  printed  antenna  array. 

2.  No  radiation  from  the  feed  network  can  interfere  with  the  main  radiation  pattern,  since  a 
ground  plane  separates  the  two  mechanisms. 

3.  No  direct  connection  is  made  to  the  antenna  elements,  so  problems  such  as  large  probe  self 
reactances  or  wide  microstrip  line  (relative  to  patch  size),  which  are  critical  at  millimeter-wave 
frequencies,  are  avoided. [6] 


Pozar  designed  a  prototype  antenna  using  small  hole  coupling  theory  and  a  cavity  model 
for  the  patch.  No  specific  bandwidth  performance  results  were  cited,  but  he  did  report  that  the 
patterns  in  the  principle  radiation  planes  were  basically  the  same  as  for  a  microstrip  antenna  with 
same  geometry.  Also  radiation  from  the  feed  side  of  the  antenna  was  negligible. 
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Sullivan  and  Schaubert  [7]  took  Pozar’s  concept  and  developed  an  exact  mathematical  model 
for  an  aperture  fed  microstrip  antenna  with  one  radiating  antenna  patch.  They  formulated  coupled 
integral  equations  by  using  the  Green’s  functions  for  grounded  dielectric  slabs.  Their  analysis 
included  all  the  coupling,  radiation  and  surface  wave  effects  for  both  substrates.  In  their  analysis 
they  invoked  the  equivalence  principle,  closed  off  the  aperture  and  replaced  it  with  magnetic  surface 
currents  Ms  located  just  above  and  below  the  ground  plane  (see  Figure  4).  To  ensure  continuity  of 
the  tangential  electric  field  across  the  aperture,  the  magnetic  currents  are  of  equal  magnitude  but 
oppositely  directed.  The  space  below  the  ground  plane  {z  <  0)  is  designated  as  region  a  and  the 
region  above  (r  >  0)  as  region  6.  The  known  incident  current  distribution  on  the  feedline  is 
the  scattered  current  distribution  on  the  feedline  is  Jj  and  the  current  distribution  on  the  patch  is 
Jp.  The  electric  and  magnetic  fields  in  each  region  can  then  be  written  as  the  summation  of  the 
contributions  from  each  source: 

=  Ea(Ji„c)  +  Ea(Jf)  +  Ea(Ms) 

=  Ha(Ji„c)  +  Ha(Jf)  +  Ha(Ms) 

E‘°t  =  Eb(Jp)-B,,(Ms) 

Htot  =  Hb(Jp)-Hb(Ms) 
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By  enforcing  the  boundary  conditions  they  obtained  three  coupled  integral  equations  for  the  un¬ 
known  currents  Jf ,  Ms  and  Jp.  The  equations  were  then  solved  using  a  Galerkin  moment  method 
procedure.  The  analysis  was  simplified  by  assuming  the  electric  currents  on  the  antenna  patch  and 
stripline  were  only  in  the  y-direction  and  the  aperture  was  electrically  shorted  and  replaced  by 
magnetic  currents  confined  to  the  x-direction.  The  expression  obtained  for  Ms  was  not  quite  exact 
because  one  material  parameter  must  be  determined  from  empirical  data.  The  authors  verified 
their  results  by  comparing  calculated  to  measured  input  impedance  for  several  different  combina¬ 
tions  of  dielectric  constant,  dielectric  thickness  and  aperture  position  and  obtained  good  agreement 
between  predicted  and  measured  performance. 

Tsao  ti  al  [8]  reported  the  results  from  the  design  and  testing  of  several  aperture  fed  stacked- 
patch  microstrip  antennas.  A  test  model  with  a  single  input  feed  network  achieved  a  19.2  percent 
bandwidth  at  3.9  GHz  with  an  input  voltage  standing  wave  ratio  less  than  two.  A  two-input-port 
model  designed-for  dual  polarization  achieved  a  23  percent  bandwidth  at  3.74  GHz  with  an  input 
VSWR  less  than  two.  No  exact  analysis  was  presented  but  the  authors  did  demonstrate  that 
an  aperture  fed  stacked-patch  microstrip  antenna  can  yield  higher  bandwidth  than  a  single  patch 
microstrip  antenna. 

Das  [9]  studied  the  general  configuration  of  printed  antennas  fed  by  planar  transmission 
lines  in  a  multilayered  substrate  geometry  by  analysis  using  a  full  wave  spectral  domain  moment 
method.  In  his  model  he  developed  a  complete  set  of  full  wave  spectral  domain  Green’s  functions 
for  a  multilayered  geometry.  One  particular  configuration  he  studied  was  an  aperture  fed  stacked- 
patch  antenna.  He  verified  his  analysis  by  using  a  waveguide  simulator  and  reported  a  bandwidth 
of  11  to  almost  17  percent  depending  upon  the  primary  patch  size. 

In  a  recent  article  Pozar  [10]  modeled  an  infinite  phased  array  of  aperture  coupled  microstrip 
patches.  His  solution  used  the  exact  Green’s  functions  for  the  dielectric  substrates  in  a  spectral 
domain  moment  method  approach.  The  coupling  of  one  aperture  in  the  infinite  array  to  its  ar>- 
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sociated  feedline  is  found,  and  an  equivalent  series  impedance  is  determined.  The  model  assumes 
no  coupling  between  the  feed  lines  themselves  as  long  as  they  do  not  come  within  a  few  substrate 
thicknesses  of  each  other.  A  waveguide  simulator  was  used  to  verify  the  theory.  The  agreement 
between  measured  and  calculated  impedance  results  was  good. 

S,S  Analysis  of  Microsirip  Antennas 

The  electromagnetic  fields  radiated  by  a  microstrip  antenna  are  found  by  solving  the  integral 
equations  for  the  currents  on  the  feed  elements  and  the  antenna  patch(s).  Mosig  [11]  describes  a 
mixed  potential  integral  equation  (MPIE)  technique  applied  to  microstrip  structures  which  can  be 
used  to  solve  for  the  currents.  In  this  method  he  uses  Green’s  functions  associated  with  the  scaler 
and  vector  potentials  which  are  calculated  by  using  stratified  media  theory  and  are  expressed  as 
Sommerfeld  integrals.  The  author  solves  the  MPIE  in  the  space  domain  to  keep  a  good  physical 
insight  to  the  problem.  The  MPIE  is  numerically  stable  and  can  be  solved  with  efficient  algorithms. 
The  formulation  of  these  equations  provides  a  powerful  and  flexible  technique  because  contributions 
from  coupling,  dispersion,  radiation  losses  and  surface  waves  are  included.  The  technique  can  be 
extended  to  multilayered  substrates  and  multiple  conductors  at  different  levels  (stacked  patches) 
by  making  suitable  modifications  to  the  Green’s  functions  and  increasing  the  number  of  unknowns. 

Additional  insight  into  using  the  MPIE  is  provided  by  an  earlier  paper  by  Mosig  and  Gar- 
diol  [12].  The  authors  develop  in  detail  the  spatial  Green’s  functions  associated  with  microstrip 
structures.  These  Green’s  functions  correspond  to  the  fields  and  potentials  created  by  a  horizontal 
electric  dipole  placed  on  the  air/dielectric  interface  and  are  expressed  as  Sommerfeld  integrals.  The 
situation  where  both  the  source  and  observer  are  located  on  the  same  plane  is  studied  extensively. 
Numerical  methods  to  evaluate  the  Sommerfeld  integrals  in  this  situation  are  reviewed,  and  several 
techniques  introduced  and  discussed  in  detail. 

Analytical  techniques  applicable  to  stratified  media  were  needed  to  develop  the  model  for  the 
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antenna  studied  in  this  thesis.  Mosig  indicated  the  MPIE  can  be  used  for  stratified  media  but  he 
derived  the  Green’s  functions  using  only  the  horizontal  electric  dipole  (HED)  and  a  point  charge 
and  did  not  derive  the  Green’s  functions  for  stratified  media  [11];  Kong  developed  the  integral 
expressions  for  the  electric  and  magnetic  fields  for  both  the  HED  and  the  horizontal  magnetic 
dipole  (HMD)  in  both  infinite  and  semi-infinite  stratified  media  [13].  As  discussed  later,  Nazar 
extracted  the  HED  and  HMD  Green’s  functions  for  stratified  media  from  Kong’s  integral  field 
expressions  [3].  The  Green’s  functions  for  the  HMD  were  needed  to  determine  the  fields  radiated 
by  the  magnetic  currents  used  to  close  off  the  aperture  in  the  ground  plane. 

Das  and  Pozar  derived  a  generalized,  two  dimensional,  spectral-domain  Green’s  function 
.completely  defining  the  field  inside  a  multilayered  dielectric  structure  due  to  a  current  element 
located  between  any  two  layers.  They  derive  their  solution  by  solving  a  “standard”  form  containing 
the  current  element  located  on  the  interface  between  any  two  layers  and  using  an  iterative  algorithrn 
to  take  care  of  additional  layers.  Another  iterative  algorithm  is  then  used  to  find  the  field  in  any 
layer  in  terms  of  the  field  expressions  in  the  two  layers  of  the  “standard”  form.  The  locations  of 
the  poles  of  the  Green’s  function  are  predicted,  and  an  asymptotic  form  is  derived  along  with  the 
asymptotic  limit.  From  their  analysis  they  draw  the  following  conclusions: 

1.  The  convergence  of  the  numerical  solution  of  the  Green’s  functions  is  much  faster  when  the 
source  and  observer  are  on  different  planes  than  when  they  are  coplanar. 

2.  The  numerical  convergence  to  the  asymptotic  limit  is  much  faster  for  thicker  dielectric  layers. 

3.  The  convergence  of  the  numerical  solution  is  much  faster  for  structures  with  a  ground  layer 
than  for  structures  without  a  ground  plane  [14]. 

Most  of  the  exact  analysis  of  microstrip  antennas  has  concentrated  on  deriving  the  Green’s 
functions  using  only  horizontal  dipole  sources  such  as  the  HED  and  HMD.  Hall  and  Mosig  [10] 
applied  the  MPIE  approach  to  vertical  monopoles  embedded  in  a  dielectric  substrate  media.  The 
vector  and  scalar  potentials  were  written  in  terms  Sommerfeld  integrals  and  evaluated  using  tech¬ 
niques  similar  to  those  developed  for  horizontal  dipoles  in  microstrip  antennas.  The  analysis  was 
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verified  by  comparing  measured  and  calculate  impedance  data  for  a  monopole  embedded  in  Teflon. 
Good  results  were  obtained  indicating  this  method  may  be  useful  for  modeling  the  coaxial  feed  for 
a  microstrip  patch  antenna. 

Barlatey  ti  al.  analyzed  a  stacked  patch  microstrip  antenna  using  the  MPIE  approach.  They 
provided  a  complete  and  rigorous  treatment  of  the  Green’s  functions  for  the  potentials  in  the  space 
domain,  created  by  a  HED  located  at  any  of  the  dielectric  interfaces  of  a  multilayered  structure.  The 
surface  waves  which  can  arise  in  such  a  structure  were  characterized  numerically.  An  approximate 
formula  for  the  value  of  the  propagation  constant  of  the  dominant  TM©  mode  was  given,  and 
conditions  necessary  to  ensure  the  absence  of  higher  order  modes  were  specified.  The  MPIE  was 
solved  with  a  Galerkin  method  of  moments  procedure.  They  show  that  the  introduction  of  auxiliary 
variables  greatly  reduces  computation  time.  Theoretical  results  for  the  input  impedance  of  a  two- 
layer  patch  agreed  well  with  measurements.  The  authors  estimate  one  configuration  would  have  a 
bandwidth  greater  than  15  percent  for  a  VSWR  of  two  or  less. 

2.3  Nazar’s  Thesis[3] 

Nazar  developed  a  theoretical  model  for  the  analysis  of  an  aperture  fed  stacked-patch  mi¬ 
crostrip  antenna.  The  mixed  potential  integral  equation  (MPIE)  approach  was  used.  The  aperture 
was  closed  by  using  opposing  magnetic  currents  on  each  side  of  the  ground  plane.  As  noted  before, 
Nazar  derived  the  stratified  media  Green’s  functions  using  both  the  HED  and  HMD.  The  necessary 
Green’s  functions  associated  wjth  the  vector  and  scalar  potentials  were  evaluated  in  the  spacial 
domain  to  keep  a  good  phj  sica/  insight  to  the  problem.  The  Green’s  functions  are  expressed  as 
Sommerfeld  integrals.  A  r.isthoc!  of  moments  technique  to  solve  for  the  currents  of  the  antenna 
is  outlined  but  no  actual  solutions  are  calculated.  The  Sommerfeld  integral  in  the  Green’s  func¬ 
tions  were  analyzed  and  found  to  have  the  following  characteristics:  complex,oscillatory  integrands; 
singularities;  surface  waves;  and  semi-infinite  integration  intervals.  He  also  developed  several  inte- 
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gration  techniques  to  deal  with  these  characteristics.  Example  calculations  for  several  of  the  Green’s 
functions  were  performed.  The  cut-off  frequencies  for  the  surface  waves  were  determined  and  he 
found  that  for  the  proper  choice  of  material  parameters  only  one  surface  wave  mode  propagates. 
Example  Green’s  functions  were  evaluated  with  sample  results  reported. 

2,4  Summary 

Microstrip  antenna  research  continually  produces  innovative  designs  and  more  elaborate  anal¬ 
ysis  techniques.  Work  by.  Pozar,  Sullivan  and  Schaubert,  Tsao  et  al,  Das,  and  Barlatey  ei  ai  have 
shown  the  feasibility  of  using  aperture  coupling  and  stacked  patches  to  increase  the  bandwidth  of 
microstrip  antennas.  Sullivan  and  Schaubert,  Mosig  and  Gardiol,  Kong,  and  Das  and  Pozar  devel¬ 
oped  analytical  techniques  to  rigorously  calculate  the  solutions  for  microstrip  antennas  in  stratified 
media.  Hall  and  Mosig  applied  the  MPIE  approach  to  vertical  monopoles  embedded  in  a  dielectric 
substrate.  Nazar  developed  the  Green’s  functions  for  an  aperture  fed  stacked-patch  microstrip 
antenna.  Solving  these  Green’s  functions  numerically  in  an  efficient  manner  will  allow  their  use  in 
a  method  of  moments  analysis  to  determine  the  antenna  currents.  The  complete  analytical  model 
for  this  antenna  can  then  be  used  to  define  frequency,  bandwidth,  and  radiation  characteristics  and 
thus  develop  criteria  for  the  efficient  design  of  this  antenna. 
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III.  Theory 


S.l  Analysis  Overview 

The  goal  of  this  analysis  is  to  find  the  induced  electric  currents  on  the  two  antenna  patches 
given  an  assumed  magnetic  surface  current  over  the  aperture.  These  currents  can  be  used  to  de¬ 
termine  the  antenna  operating  characteristics  such  as  resonant  frequency,  input  impedance,  band¬ 
width,  and  radiation  pattern.  The  various  antenna  currents  are  illustrated  in  Figure  5.  The  electric 
currents  are  labeled  as  for  the  feedline  current;  Jf  is  the  feedline  scattered  current;  J2  is  the 
primary  patch  current;  and  J3  is  the  parasitic  patch  current.  The  equivalence  principle  [15]  can  be 
used  to  close  off  the  aperture  and  replace  it  with  opposing  magnetic  surface  currents  Mj  located 
just  above  and  below  the  ground  plane.  The  opposing  magnetic  currents  ensure  continuity  of  the 
tangential  electric  field  through  the  aperture  [7].  Though  only  x-directed  electric  and  y-directed 
magnetic  currents  are  shown  in  Figure  5,  the  analysis  includes  all  currents  in  the  xy-plane.  The  total 
electric  and  magnetic  fields  in  regions  a  and  6  can  be  written  as  the  summation  of  the  contributions 
from  each  source: 


=  Ea(Jinc)  +  Ea(Jf)  +  Ba(Mi) 

(1) 

=  Ha(Ji„c)  +  Ha(Jf)  +  Ha(Mi) 

(2) 

=  Eb(J2)  +  Eb(J3)-Eb(Mi) 

(3) 

4°*  =  Hb(J2)  +  Hb(J3)-Hb(M8) 

(4) 

In  this  analysis  a  somewhat  simplified  model  is  used.  Rather  than  using  as  the  known 
input,  the  magnetic  surface  current  across  the  aperture  is  assumed  known.  Thus  only  the  two 
patch  currents  J2  and  J3  need  be  solved  for.  To  determine  these  currents,  only  expressions  for  the 
electric  field  need  to  be  derived. 
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The  electric  fields  in  the  problem  can  be  expressed  using  the  mixed  potential  integral  equation 
(MPIE)  and  solved  in  the  space  domain.  Harrington  discussed  using  both  vector  and  scalar  poten¬ 
tials  (mixed  potentials)  to  solve  scattering  and  antenna  problems  [16].  Mosig  applied  the  MPIE  to 
microstrip  antennas  [11]  and  Nazar  applied  them  to  the  aperture  fed  stacked  patch  structure  [3]. 

The  electric  fields  due  to  the  electric  and  magnetic  sources,  respectively,  are  derived  from  the 
scalar  and  vector  potentials 


E  =  -jwA-VV  (5) 

E  =  --VxF  (6) 

£ 

where  A  is  the  magnetic  vector  potential,  V  is  the  electric  scalar  potential,  F  is  the  electric  vector 
potential,  and 

Equation  (5)  assumes  only  electric  sources  are  present,  and  Eq.-(6)  assumes  only  magnetic 
sources.  The  vector  and  scaler  potentials  can  be  expressed  using  the  corresponding  Green’s  function 
as  superposition  integrals  of  the  current  and  charge  densities. 


A(p)  = 

jfGA(plp')-J(p')ds' 

(7) 

VM  = 

/  Gq(p!p')q(p')ds' 

(8) 

F(P)  = 

/  Gf(/>|p')  •  M(p')ds' 

Js’ 

(9) 

with  the  dot  in  Eqs  (7)  and  (9)  indicating  a  dot  product  of  the  dyadic  Green’s  function  with  the 
vector  surface  current.  The  radial  positions  of  the  observer  and  source  respectively,  from  the  z-axis 
arc  represented  by  p  and  p'.  The  continuity  equation  relates  the  electric  current  and  charge  density. 

V-J-fjwq  =  0  (10) 
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The  electric  field  in  region  b  (the  half  space  into  which  the  antenna  radiates)  is  described  by 
substituting  Eq  (7)  through  (10)  into  Eqs  (5)  and  (6).  Two  coupled  integral  equations  are  then 
obtained  for  the  two  unknown  currents  J2  and  J3  by  enforcing  the  boundary  conditions  of  Table  1. 
The  remaining  boundary  conditions  such  as  conditions  on  the  normal  and  tangential  fields  at  the 
dielectric  interfaces,  are  incorporated  in  the  construction  of  the  Green’s  functions. 

Table  1.  MPIE  Boundary  Conditions  [3] 

1)  E*""  =  0  on  antenna  patch  1 

2)  E‘“"  =  0  on  antenna  patch  2 


In  developing  a  mathematical  model  the  first  step  is  to  determine  the  expressions  for  the  vector 
and  scalar  potentials  from  which  the  Green’s  functions  are  constructed.  The  MPIEs  satisfying  the 
boundary  conditions  in  Table  1  are  then  formulated  and  can  be  solved  for  the  current  and  charge 
distribution  by  a  moment  method  solution.  Throughout  the  analysis  a  time  dependence  of  is 
assumed. 

S.2  Vector  and  Scalar  Potentials 

By  definition  the  Green’s  functions  are  potentials  created  by  a  unit  source  which  is  an  electric 
or  Hertz  dipole  located  on  of  the  patches  or  an  equivalent  magnetic  source  located  across  the 
aperture  in  the  ground  plane.  The  superposition  principle  applies  in  a  linear  system  and  the 
potentials  of  any  finite  source  can  be  determined  by  representing  the  source  as  a  continuum  of 
elementary  dipoles  and  integrating  the  contributions  of  all  the  elementary  sources  [17].  The  method 
of  Green’s  functions  for  arbitrary  microstrip  sources  was  developed  by  Mosig  and  Gardiol  [12,  17]. 
However,  they  only  constructed  Green’s  functions  for  a  single  dielectric  layer  antenna,  with  only 
electric  sources  on  the  dielectric-air  interface  and  the  dielectric  permeability  fixed  at  the  free  space 
value.  Nazar  [3]  extended  these  results  to  model  an  aperture  fed  stacked-patch  structure  with 
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two  dielectric  layers  with  different  permittivities  and  permeabilities;  with  electric  sources  at  the 
dielectric-dielectric  and  dielectric-air  interfaces;  and  with  equivalent  magnetic  sources  at  the  ground 
plane-dielectric  interfaces. 

S.S.l  HED  ai  Interface  Sb  [S],  The  first  step  in  analyzing  this  structure  is  to  determine  the 
magnetic  vector  potential  created  by  a  HED  having  a  unit  moment  {Idx  =  1  A  m)  along  the  x-axis 
at  the  interface  of  dielectrics  lb  and  2b  (see  Figure  6)  [3].  The  resulting  Green’s  functions  can  be 
used  in  Eqs  (7)  and  (8)  along  with  Eq  (5)  to  obtain  the  E  fields  of  the  structure  for  an  arbitrary 
distribution  of  sources  on  interface  2b.  The  x  and  z  components  of  the  magnetic  vector  potential 


Figure  6.  HED  at  Interface  2b  (dielectric  2b  is  transparent  for  clarity)  [3]. 
in  each  dielectric  and  free  space  can  be  expressed  in  the  form  of  Sommerfeld  integrals  [17] 


Ax12(p) 
■^21 2  (p) 


j  Ho^\Ap)a52sinh(uibz)dA 
cos  <{)  j  H^^^(Ap)a^2  cosh(uibz)dA 


(11) 

(12) 
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Ax22(p)  =  [b^2S»nh(u2bz)  +  Cx2  sinh(u2bz)]  dA  (13) 

Az22(p)  =  cos  ^  y  H^i^^(Ap)  [b^2  sinh(u2bz)  +  c^2  sinh(u2bz)]  dA  (14) 

Ax32(p)  =  yH^*\Ap)d^2exp(-U3bz)dA  (15) 

Az32(/>)  =  cos^yH^j*)(Ap)d^2exp(-U3bz)dA  (16) 


where  ^21  ^21  ^x2i  ®^c.  are  unknown  coefficients  to  be  solved  for,  p  is  the  radial  distance  |p|  between 
the  source  HED  and  the  observer  point  and 

uib  =  \/A2-kfb,  U2b  =  ^A2-k^b>  “3b  =  ^AZ-k^b 

with  ku,  k2b,  and  kai  being  the  wave  numbers  in  dielectrics  lb,  2b  ,  and  3b,  respectively.  Although 
kzh  is  assumed  to  be  the  free  space  value,  ko,  it  is  referred  to  as  jt3i  for  notation  consistency.  The 
wave  numbers  are  defined  as 

kib  =  uy/enfin,  k2b  =  k3b  =  Uy/eatfizb 

with  eu,2i,36  and  Pi6,2»,3»  being  the  permittivities  and  permeabilities  of  the  corresponding  medium. 
As  an  example  of  the  notation,  Axj2  is  defined  as  the  x-directed  magnetic  vector  potential  in  region 
6,  dielectric  lb  (0  <  z  <  bib)  for  a  HED  at  interface  2b.  The  remaining  vector  potentials,  A^,*  nm 
are  defined  similarly.  A  definition  of  the  integration  path  C  along  with  some  background  of  the 
development  of  Eqs  (11)  -  (16)  can  be  found  in  the  works  by  Nazar  [3],  and  Mosig  and  Gardiol  [17]. 

Nazar  developed  the  necessary  boundary  condition  to  solve  for  the  coefficients  of  Eqs  (11)  - 
(16).  The  relevant  results  applicable  to  this  work  are  summarized  as  follows: 

@  z  =  bib 


Aj12  —  Ax22 


(17) 
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(18) 


1  _ 1  3Ay22  _  ^(P) 

fin  dz  .  fiih  dz  iirp 
■^zl2  _  Az22 


Pl6  P2l 

eihPih  dz 


1  M^22 

£24/121  dz 


=  COS(^ 


\e26P2i 


_1_'\  ^Ax12 
eiiPi6/  3p 


(19) 

(20) 


@  Z  =  b2b 


■^x22  ~  ■^x32 

1  3A^22  _  1  ^Ax32 
P26  dz  flat,  dz 
A722  _  -^232 
P2i  P36 

1  aA^22  1  ^A^32  .  f  1 

£26/124  02  eahpab  dz  \c3hP3h 


^  \  ^■^x32 

E2I.P2i/  dp 


(21) 

(22) 

(23) 

(24) 


along  with  the  Sommerfeld  radiation  condition 


lim  r  =  ( 

r-»oo 


dA 


+  jkA  1=0 


(25) 


where  A  satisfies  the  homogenous  Helmholtz  equation 


(V^  +  k=^)  A  =  0 


(26) 


The  branch  interpretation  [18]  of  Usb  =  -  kg^  is  determined  by  the  Sommerfeld  radiation 

condition  as 


U3b  = 


U3b  =  j  \/A2  -  k|b  , 


|A1  >  k3b 

1A|  <  ksb 
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The  term  5(p)/(27rp)  in  equation  (18)  represents  a  unit  current  source  {Idx  =  1  A  m)  on  interface 
lb,  where  S{p)  is  the  Dirac  delta  function  which  can  be  expressed  as  [17] 


(27) 


Using  the  boundary  conditions  defined  and  Eqs  (11)  -  (16),  eight  equations  are  formed  and  solved 
for  the  unknown  coefficients  [3]  yielding 


Ax12(p) 

■^x22(/’) 

Az32(p) 


=  J_  [e^^\Xp)  [NL2Wsinh(u2bz)  +  N^,2(A)sinh(u2bz)1  ^ 

4it  Jf.  L  Dt('^) 

_ »  TxtK  /\\_»  V/..  _\  I  xtK  /\\_»  1./..  -.\1 


COS(j) 

47r 


^Hf(Ap) 


N^^2(A)sinh(u2bz)  +  N^2(A)sinh(u2bz) 


D.^(^)D^(A) 


(28) 

(29) 

(30) 

(31) 

(32) 

(33) 


with  the  N’s  representing  the  numerators  of  5^2i^2>^x2i^22>  represent  the  denomina¬ 

tors.  The  zeros  of 


,  [pH2U3b  +  UlbCOth(bibUlb)]U2bCOsh(U2b(b2b-blb))  ,  ^ 

D^(A)  =  (34) 

+  [pn2U2b  +  PW3Uib«3bCot(bibUib)]  sinh(u2b(b2b  -  bib)) 

.  [en3U3b  +  uib  tanh(bibUib)]  U2b  cosh(u2b(b2b  -  bib)) 

D^(A)  =  (35) 

+  [e6i2U^b  +  «l23UibU3b  tan(bibUib)]  sinh(u2b(b2b  -  bib)) 


define  the  surface  wave  poles  in  the  dielectrics  where 


The  frequencies  at  which  the  TE  and  TM  surface  waves  propagate  are  deterihined  by  the  roots 
of  De(A)  and  Dm(A)  [3].  The  surface  waves  represent  power  propagating  along  the  surface  of  the 
dielectric  layers,  instead  of  radiating  into  space.  De(A)  and  D^(A)  are  related  to  the  reflection 
factors  on  a  microstrip  structure  for  an  incident  perpendicular  polarized  wave  (TE)  and  parallel 
polarized  wave  (TM),  respectively  [17].  The  N  terms  are  given  in  Appendix  A. 


S.S.S  HED  at  Interface  Sb  [Sj.  The  case  for  a  HED  located  on  the  dielectric  2b  and  dielectric 


3b  (free  space)  is  shown  in  Figure  7. 


Figure  7.  HED  at  Interface  3b  [3]. 


The  X  and  z  components  of  the  magnetic  vector  potential  in  each  dielectric  and  free  space  [3] 


can  be  expressed  as 


AxisCp)  =  j^H[,*^(Ap)a^3sinh(uibz)dA  (36) 

A*i3(p)  =  cos^  yHf^(Ap)a^3Cosh(uibz)dA  (37) 

Ax23(p)  =  j^H[,*^(Ap)  [b$3sinh(u2bz)  +  Cx3  sinh(u2bz)]  dA  (38) 

AxJ3(p)  =  cos^yHf^(Ap)  (b^3sinh(u2bz)  +  c^sinh(u2bz)]  dA  (39) 

Ax33(p)  =  j^H^^^(Ap)d^3exp(-U3bz)dA  (40) 

A233(p)  =  yHf^(Ap)d^3exp(-U3bz)dA  (41) 


The  boundary  conditions  needed  to  solve  for  the  coefficients  of  equations  (36)  -  (41)  are: 
@  Z  :=  bib 

Ax13  =  Ax23 
^  ^Ax13  _  1  ^Ay;3 

Pi  6  dz  P2t  dz 
^213  _  Az23 
Pll  P26 

_L_£^  _  =  C0S6  (— _ —]  ^ 

cnnib  dz  e2hH2b  dz  \e2^P2^  euPU/  dp 


@  z  =  b2b 


(42) 

(43) 

(44) 

(45) 


Ax23  —  Ax33 

1  9A^23 _ 1  dA^23  _  ^(p) 

P2b  dz  pzb  dz  2np 

■^223  _  -^233 

P2l  UZb 

_  _J_d^  ^  _ 1_\  aA^33 

e26P26  dz  ezbUzb  dz  ^  \e36P36  ^2bli2b)  dp 


(46) 

(47) 

(48) 

(49) 
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along  with  the  Sommerfeld  radiation  condition  of  Eq  (25).  As  before,  the  term  6{p)/2itp  in  Eq  (47) 
represents  a  unit  current  source,  but  this  time  on  interface  2b.  Solving  for  the  unknown  coefficients 
r^ults  in  the  following 


AxIsCp) 

A*i3(p) 

Ax23(p) 

A223(p) 

Ax33(p) 

A233(p) 


J_  jfH(2)(Ap)  ^N^x3(^)sinh(u2bz)  +  N^„3(A)sinh(u2bz)j 

^  r^b73Wsinh(u2bz)  +  N^23(A)sinh(u2bz)'| 

’T  Jc  *  '  L  D.^(A)D  m  W 


cos 

4ff 


477 
COS^ 

4^ 


Jnf\xp) 


Nb23(^)sinh(u2bz)  +  N^23(A)sinh(u2bz)  j 
-**’‘^--)exp(-U3bz)dA 

NL(A) 


D."(A)D^(A) 


exp(-U3bz)dA 


(50) 

(51) 

(52) 

(53) 

(54) 

(55) 


The  N  terms  of  equations  (50)  -  (55)  are  included  in  Appendix  A. 


S.2.3  HMD  at  Interface  lb  [3],  The  electric  fields  in  region  b  due  to  the  equivalent  magnetic 
current  on  the  ground  plane  (see  Figure  8)  are  found  with  Eq  (9)  along  with  Eq  (6).  The  first  step 
in  the  analysis  is  to  find  the  fields  created  by  a  HMD  along  the  z-axis  at  the  ground  plane  with 
a  unit  moment  Vdz  =  1  V  m.  Using  the  same  procedure  as  for  the  electric  sources,  the  z  and  z 
components  of  the  electric  vector  potential  in  each  dielectric  and  free  space  of  region  6  are  expressed 
as 


Fxii(p)  =  jfl4^^(A/7)[a$,sinh(uibz)  +  exi  cosh(uibz)]  dA  (56) 

F^iCp)  =  cos()ij^Hi^^(Ap)a^^isinh(uibz)dA  (57) 

Fx2i(p)  =  ^  H^^^Ap)  [bxi  sinh(U2b2)  +  Cxi  cosh(u2bZ)]  dA  (58) 

Fz2i(p)  =  cos^yHf^(Ap)  [bjisinh(u2bz)  +  C2iCosh(u2bz)]  dA  (59) 

Fx3i(p)  =  j^H^^^(Ap)d^iexp(-U3bz)dA  (60) 
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z 


Figure  8.  HMD  at  Interface  lb  (dielectrics  lb  k  2b  are  transparent  for  clarity)  [3]. 

F23i(/»)  =  cos^j^Hi^^(Ap)d^iexp(-U3bz)dA  (61) 


The  boundary  conditions  needed  to  solve  for  the  unknown  coefficients 
@  z  =  0 

1  dFtn  ^  6{p) 

cij  dz  2np 


are: 


(62) 


@  z  =  bib 


F^u=F,V 


£lJ  dz  C2J  dz 

!k  =  !k 

eit  C2i 

^  _ ^ 

eu/iij  dz  e2ii^2i  oz  \£2ili2b  snpnj  dp 


(63) 

(64) 

(65) 

(66) 
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®  z  =  b2b 


—  Fx31 


(67) 


1  aF^ai  ^  1 
e2i  dz  C31  dz 

!k  =  !k 

C2i 

e3i/<3i 


=  cos^ 


\C3hlJ^3i 


_J_)^ 

£21/^21/  5p 


(68) 

(69) 

(70) 


along  with  the  Sommerfeld  radiation  condition,  Eq  (25),  with  A  replaced  with  F,  As  before,  the 
term  -S{p)/{2irp)  in  Eq  (62)  represents  the  unit  magnetic  current  source  on  the  ground  plane 
(interface  lb).  The  solutions  of  Eqs  (56)-(61)  are  as  follows 


F^n(p) 

F^^^(P) 

F^2i(p) 

F*"2i(p) 

F^3i(p) 

Fz^i(p) 


1  ^  \  euA  .  N^,i(A)cosh(uibz)‘ 


cos<j> 

4jr 

47r 


j[HW(Ap) 


NLi(A) 


•  sinh(uibz)dA 


De"{A)D^(A) 

[NLi(A)  sinh(u2bz)  +  N^xi(A)  cosh(u2bz)] 
D^(A) 


dA 


dA 


cos<;i  /‘„(2),^  rN^^,(A)sinh(u2bz)  +  N^zi(A)cosh(u2bz)l 

Jc  [  dF(a®a) 


45r 
cos^ 
4;r 


NL(A). 
D^{A) 
NSz»(A) 


^  h[,^\Ap)  ^g»^^exp(-U3bz)dA 


De"(A)D^(A) 


exp(-U3bz)dA 


(71) 

(72) 

(73) 

(74) 

(75) 

(76) 


The  numerator  terms  are  listed  in  Appendix  A. 


S.S.4  The  Eleciric  Scalar  Potential  [3].  To  determine  the  electric  field  of  Eq  (5)  the  electric 
scalar  potentials  are  needed.  The  potentials  can  be  found  using  the  Lorentz  gauge  condition  for 
electric  sources: 


V  •  A  +  jwepV  =  0 


(77) 


When  applied  to  the  magnetic  vector  potentials  previously  calculated,  Eq  (77)  will  yield  the  po¬ 
tentials  for  two  opposite  electric  point  charges.  To  determine  the  Green’s  functions  of  Eq  (8),  the 
potentials  of  a  single  point  charge  are  needed.  The  electric  potential,  Vj,,  of  a  single  point  charge 
can  be  found  using  the  electrostatic  relationship  linking  an  electric  point  charge  and  a  dipole  po¬ 
tential  [17]: 

v  =  -^  (78) 

V,  can  be  found  within  an  arbitrary  constant  by  integrating  V  over  x.  The  value  of  the  constant  is 
unimportant  since  the  gradient  of  is  taken  to  calculate  the  electric  field. 


S.S.4‘1  Scaler  Potential  for  HED  at  Interface  Sb  [Sj.  Applying  Eq  (77)  to  Eqs  (28) 


and  (29)  yields 


^Ax12 

9x 


dz 


(79) 


resulting  in 


V^2(P)  = 


cos^  f  „(2)  [a  _  N^^2(A) 

ATriuenmJc  ‘  D^(A)  ‘*'D^(A)D^(A) 


sinh(uibz)dA 


(80) 


As  previously  stated,  this  potential  is  for  two  point  charges.  Applying  Eq  (78)  to  Eq  (80)  results 


in 


''hM  [w- 

The  potentials  in  the  other  mediums  can  be  found  similarly: 


-r-  ■  ■  _  “at ...  _ , 

"  D.^(A)Db(A)J 


I’^cxaC-^)  _  Hat 
tWiT  a  nb/t-vribr 


Dp)  ^  De^(A)DS,(A)J 


sinh(«2iz) 

cosh(u262) 


dA(82) 


yb  („\ _ ^  \  f^^dxaW  ^3b  N^22(A) 

-  4nju>e3mi  Jc  °  ^  D,^(A)  A  Db(A)D^.(A) 


exp(-U3bz)dA  (83) 
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S.2.4.S  Scalar  Poieniials  for  HED  ai  Interface  Sb. 


V^13(P)  = 


V$23(P)  = 


V$33(P)  = 


1 


AirjucihtJiit, 


(84) 


1 


Avjucihm 


1 


Airjuezifizi, 


[^)  -  5|^]  -h(u.,.)aA 

i  -  TsU^a)] 


S.S  Green’s  Functions  Construction 

The  electric  fields  in  any  medium  may  be  found  using  Eqs  (5)  and  (6)  along  with  the  vector 
potentials  and  the  electric  scalar  potentials.  However,  only  the  tangential  eclectic  fields  listed  in 
Table  2  are  needed  for  this  analysis. 


Table  2.  Necessary  Tangential  Fields  [3] 


HED  at  interface  2b  (Patch  1): 

HMD  at  interface  lb  (Aperture): 

1)  E*®"  at  interface  2 

1)  E*®"  at  interface  2 

2)  E*®"  at  interface  3 

2)  E‘®"  at  interface  3 

HED  at  interface  3b  (Patch  2): 

1)  E*®"  at  interface  2 

2)  E‘®"  at  interface  3 

S.SJ  E  Fields  Due  to  Electric  Sources  [Sj.  Equation  (5)  is  used  to  calculate  the  E  fields 
due  to  electric  sources.  The  Green’s  functions  for  the  case  of  an  electric  source  on  interface  3b  and 
an  observer  on  interface  2b  are  used  to  demonstrate  the  form  of  the  equations. 
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The  tangential  E  field  at  interface  2b  is  found  by  putting  equations  (52)  and  (85)  into  equation 
(5),  yielding 

=  [-jw  (A^^x  +  A%y)  -  (87) 

where  Ay 23  is  the  magnetic  vector  potential  for  a  y-directed  current  source  and  is  equal  to  Eq  (52). 
Since  only  a  unit  point  source  located  a  p  =  0  is  assumed  initially,  the  integrals  of  Eqs  (7)  and  (8) 
over  the  point  source  results  in 


G^&(o)  =  A,^3(p) 

(88) 

G^^/3(p)  =  Ab23(p) 

(89) 

G^^siP)  =  V,^23(p) 

(90) 

where  GassCp)  is  defined  as  the  magnetic  vector  potential  in  the  a?-direction  (first  x  in  the  super¬ 
script)  on  interface  2b  for  an  infinitesimal  x-directed  current  source  (second  x  in  the  superscript) 
on  interface  3b,  with  an  observer  at  radial  position  p.  The  Green’s  functions  in  Eqs  (89)  and  (90) 
are  similarly  defined.  By  taking  the  transverse  gradient  in  cylindrical  coordinates  and  using  the 
identity  p  =  xcos ^  +  y sin <{>,  Eq  (87)  becomes 


E‘2“3"(p)  = 


•  ^  ^®q23(p)  j 

-JwGX'aaCp) - —  cos  <f> 


x  + 


-j^G^)  - 


y  (91) 


The  Green’s  functions  only  depend  on  the  radial  separation  p  =  |p|  between  the  source  and  observer 
and  not  the  relative  angular  position.  Equation  (91)  gives  the  expression  for  the  tangential  E  field 
anywhere  on  interface  2b  for  a  point  source  at  p  =  0  on  interface  3b.  For  generalized  source  and 
observer  points  (p  0),  p  and  4>  become  R  and  C  (see  Figure  9) 


R=  |p-p'| 
^  =  sin-' 


(92) 
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Then  for  the  generalized  case  for  a  distribution  of  sources  over  interface  3b,  Eq  (91)  becomes 


Figure  9.  Geometrical  Relationships  for  x-Directed  HED  at  Arbitrary  p'  [3,  17]. 


By  properly  interchanging  the  Green’s  functions  and  source  terms  in  Eq  (93),  all  the  tangential  E 
fields  produced  by  the  HEDs  listed  in  Table  2  can  be  found.  It  appears  that  Nazar  inadvertently 
left  out  a  factor  of  jw  from  the  denominator  terms  of  Gq22(R))Gq23(R),Gq33(R),  and  Gq32(R).  A 
complete  list  of  all  Green’s  functions  can  be  found  in  Appendix  B. 


3.S.2  E  Fields  Due  io  Magneiic  Sources  [Sj.  The  tangential  E  fields  at  the  dielectric  in¬ 
terfaces  due  to  the  magnetic  source  on  the  ground  plane  are  found  with  Eq  (6).  An  example  of 
calculating  one  of  those  fields  will  be  presented. 
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The  tangential  E  field  at  interface  3b  due  to  an  infinitesimal,  a:-directed  magnetic  current 
source  on  interface  lb  (ground  plane)  can  be  found  by  expressing  F31  in  cylindrical  coordinates 

Fp3i(p)  =  ^  exp(-U3bz)dA  (94) 

(95) 

(96) 


and  then  applying  Eq  (&)  in  cylindrical  coordinates  results  in 

(97) 

(z=b,i,) 

The  X  and  y  components  of  the  E  field  can  be  determined  by  using  p  =  xcos^  +  y  sin^  and 
<})  =  -xsin^  +  y  cos^,  resulting  in: 


1 

/l5Fb3i 

£31 

Ip  H 

dz 

r+l  sz 

dp 

sin  2^ 
4ff£2l 


47re26  1 


V 


+5  /c  Ho^^(Ap)A-j^b^j(^^^  exp(-U3bb2b)dA 

-jL  Db^AjD^(A)  gxp(-M3t620dA 

-  /c  Ho*^(Ap)u3b^^^^  exp(-U3bb2b)dA 
-  cos^  ^  H^^^(Ap)A^^z^^  exp(-U3bb2b)dA 


(98) 


(99) 


where  £3^’*  designates  a  y-directed  electric  field  for  an  z-directed  magnetic  current  element.  As 
before,  the  integrals  of  Eqs  (98)  and  (99)  only  depend  on  the  radial  distance  between  the  source 
and  observer  points.  For  an  arbitrary  distribution  of  z-directed  magnetic  current  elements  let  p 
and  <j)  become  R  and  C  as  related  by  Eq  (92).  Then  Eqs  (98)  and  (99)  become 


E^r(p) 


(100) 


(101) 


E^r(p)  =  -jf^G^l(R,C)M^i(pOds' 


where  -Ge3i(R)C)  equals  the  right-hand-side  of  Eq  (98)  and  -G|3i(R,C)  equals  the  right-hand-side 
of  Eq  (99)  for  arbitrary  source  and  observer  locations. 

To  determine  the  E  fields  for  j^directed  magnetic  current  sources  replace  ^  by  ^  -  ff/2  in 
Eqs  (94)  to  (96)  [3]  and  repeat  the  same  procedure  used  to  derive  Eqs  (98)  and  (99).  The  x  and 
y  components  of  the  E  fields  at  interface  3b  due  to  an  infinitesimal,  y-directed  magnetic  current 
source  on  interface  lb  are  then 


sin  2^  J 
4jr£26  I 


-I  Sc  exp(-U3bb2b)dA  ' 

Sc  ^P^(AP)  pbf  jijb  {-X)  exp(-U3t62t)dA 
Sc  (AP)  exp(-U3bb2b)dA 

+  Sc  Ho*^(Ap)u3b^^?^  exp(-U3bb2b)dA 

+  sin^  (f>  X.  h|,’^^(Ap)A  exp(-U3bb2b)dA 


(102) 


(103) 


which  are  similar  to  Eqs  (98)  and  (99).  For  an  arbitrary  distribution  of  y-directed  magnetic  current 
elements,  the  E  fields  of  Eqs  (102)  and  (103)  are  represented  as 


^tViP)  =  -j[,G|^3"(R.C)M^i(p')ds'  (104) 

E^r(p)  =  -  (105) 

The  tangential  E  fields  can  be  written  more  compactly  as 

Etr  =  -j^^  5e3,(R,  C)  •  M?(p0ds'  (106) 
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where 


(107) 


G^KR.Oxx  G^^,iR,OW 

G|^3l(R.C)y5f  G^y3y(R,0yy 

The  tangential  E  fields  at  interface  2b  due  to  the  magnetic  surface  currents  on  interface  lb 
can  be  found  by  expressing  Fji  in  cylindrical  coordinates  and  repeating  the  procedure  for  z  =  bib 
[3]  resulting  in 

5ki(R,  C)  •  Mb(p')ds'  (108) 

The  complete  set  of  Green’s  functions  for  this  analysis  are  in  Appendix  B. 

S.4  Integral  Equations 

The  expressions  for  the  E  fields  are  integral  equations  involving  the  derived  Green’s  functions 
and  the  unknown  current  and  charge  distributions.  The  unknown  distributions  on  the  patches  are 
found  by  enforcing  the  boundary  conditions  of  Table  1.  To  comply  with  the  boundary  condition 
E‘®"  =  0  on  patch  1,  the  sum  of  the  tangential  electric  fields  due  to  the  sources  on  patch  1,  patch 
2  and  the  equivalent  sources  on  the  aperture  must  equal  zero.  Using  the  derived  Green’s  functions 
and  Eq  (5),  this  can  be  expressed  as 

-jw  /  Ga22(pIp')  •  J2(p')ds' 

Jpi 

-jw  /  GA23(p|p')-J3(p')ds' 

JP2 

+  /  Ge2i(p|p') 

JAp 

with 

Ga22(p|p0  = 


Ga22(pIp')xx  0 

0  G^y/2(^|p')yy 


(110) 


-  v‘/  G^M'Mpys' 

JPI 

Jpl 

■  Mi(p')*'  =  0  (109) 


Ge3i(RiO  — 
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Ga23(pIpO ’s  similarly  defined  and  Mi(p')  is  the  assumed  magnetic  current  across  the  aperture. 
On  patch  2 


-jw  /  Gx32(p1p')  •  h(p')is'  I  Gq32(Plp')q2(p')<Js' 

Jpi  Jpi 

-jw  /  gL3(pIp')  •  J3(p')<ls'  -  /  Gq33(p|p')q3(p')ds' 

JP2  Jpi 

+  /  ii(p|p')-Mi(p')ds'  =  0  (111) 

JAp 

The  unknown  current  and  charge  distributions  in  Eqs  (109)  and  (111)  can  be  solved  using  a  method 
of  moments  approach  which  converts  the  integral  equations  to  matrix  equations. 

3.5  Approximation  Theory 

Given  a  set  of  discrete  values  for  a  function,  approximations  for  intermediate  values  can  be 
determined  by  two  separate  but  related  processes  (19).  If  the  values  are  assumed  accurate  then 
an  interpolation  scheme  can  be  used.  This  involves  determining  a  function  whose  curve  will  pass 
through  each  data  point  or  a  chosen  subset  of  those  points.  If  the  discrete  data  are  thought  to 
contain  errors  then  a  curve  fitting  process  can  be  used  to  try  and  average  out  these  errors. 

The  Green’s  function  data  to  be  approximated  was  assumed  accurate,  so  the  first  thought 
might  be  to  use  an  interpolation  technique.  The  numerical  data  for  each  Green’s  function  consisted 
of  50  data  sets  which  would  require  an  interpolation  polynomial  of  order  49  if  all  the  data  sets  were 
used.  The  process  of  determining  the  coefficients  of  this  size  polynomial  would  only  be  trading  one 
computationally  intense  process  for  another.  Using  only  a  subset  of  the  data  would  be  possible 
if  one  knew  a  priori  which  data  sets  should  be  used.  No  previous  knowledge  of  the  data  to  be 
approximated  was  assumed,  so  ail  the  data  was  used  in  determining  the  approximations.  Thus 
to  use  all  the  data,  a  curve  fitting  method  was  chosen  to  determine  a  polynomial  or  functional 
approximation  of  hopefully  much  lower  order.  A  linear  least  squares  method  was  chosen  for  use  in 
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determining  the  coefficients  of  the  approximation  functions. 


In  the  least  squares  method  a  form  of  function  or  a  polynomial  is  assumed  and  the  coefficients 
are  solved  for  based  upon  the  goal  of  minimizing  the  square  of  the  error  between  the  original  function 
and  the  approximation.  An  outline  of  the  method  as  presented  by  Stark  [19]  follows. 

Given  a  table  of  n  + 1  points,  an  approximating  function  of  the  form 

Pm(*)  =  am9m{x)  +  +  •  •  •  +  aiffi(x)  +  aoifo(a:)  (112) 

is  assumed,  where  the  functions  (/m(®)i  •  •  •  are  some  known  functions  of  x.  The  approximating 
function  is  sought  to  match  the  n  +  1  given  points  in  such  a  way  that  if  a  set  of  deviations  were 
formed,  one  for  each  tabulated  x,-, 

50  =  Pm(®o)-/(xo) 

51  =  Pm(*l)-/(®l) 

52  =  Pmi^T)  -  f{x2) 

Sn  —  Pm(®n)  ~  f{Xn) 

or,  in  general  terms,  the  set  of  deviations 

Si  =  Pm{xi)  -  fixi)  fori  =  0,1, 2,...,  n  (113) 

then  the  sum  of  the  squares  of  these  n  +  1  deviations  should  be  a  minimum: 

n 

'^{Si)^  =  a  minimum. 

«=o 
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(114) 


The  summation  in  Eq.  (114)  can  also  be  written  as 


=  E[Pm(x.)-/(*.)]^  (115) 

t=0  t=0 

=  52[amffm(*.-)  +  "-  +  aoffo(*i-)-/(*0]^  (116) 

i=0  i=0 

The  object  is  to  minimize  the  right-hand  side  of  Eq.  (116).  Since  the  Xi  and  /(x,)  are  given  in  the 
original  table,  and  the  </,•(*)  are  known  functions  which  can  be  evaluated  at  the  given  Xj,  the  only 
unknowns  in  Eq.  (116)  are  the  coefficients  Of.  Thus  Eq.  (116)  can  be  looked  on  as  a  function  of 
the  m  -1- 1  coefficients  a,-  with  all  the  other  terms  being  known  constants. 

To  minimize  a  function  its  first  derivative  is  set  equal  to  zero  and  the  equation  is  solved  for 
the  unknown  variable.  In  this  case  the  function  to  be  minimized,  Eq  (116),  is  actually  a  function  of 
the  m-h  1  unknowns  aj.  The  first  step  in  solving  for  the  unknowns  is  to  find  the  partial  derivative 
of  the  summation  with  respect  to  each  of  the  m  -f  1  variables  Oj,  and  set  each  of  these  derivatives 
equal  to  zero: 


»=o 


0 

0 


n 


0. 


or,  in  general  terms 


g  n 

T“5Z(^‘)^  =  0  fori  =  0,1,2,...,  m. 


(117) 
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The  series  has  a  finite  number  of  terms,  and  so  the  derivative  of  a  sum  is  equal  to  the  sum  of  the 
derivatives  of  the  terms.  Thus  Eq  (117)  can  be  written  as 


=  0  fori  =  0,1, 2 . (118) 

i=0 

The  differentiation  of  any  6,-  with  respect  to  Cj  is  simple  when  noting  that 

Si  =  Pm{Xi)-f{Xi) 

=  <im9m{Xi)  +  am-lflm-l(®«)  + - h  <^j9j{^i)  +  ’  ’  ’ 

+aigi{xi)  +  ao9o{xi)  -  /(x<).  (119) 


When  differentiating  with  respect  to  aj,  every  other  term  in  Eq  (119)  is  assumed  constant  except 
the  term  containing  aj,  so  that  all  the  other  terms  drop  out  upon  differentiating,  and 

d6'  d 

— i-  =  0  +  0  +  *"+'^ — oypj(xi)  +  0  + - 1-0 

C/dj  OQj 

H  =  «(..)•  m 

Now  substituting  Eq  (120)  into  Eq  (118),  m+  1  equations  are  obtained: 


n 

X^%(®i)  =  0  fori  =  0,1,2, ,.,,m 

t=0 


(121) 
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Substituting  Eq  (113)  into  Eq  (121,  results  in 


X^bm(a:i) -/(*«•)]  =  0  for  j  =  0,1,2,...,  m,  (122) 

i=0 


and  then  substituting  Eq  (112)  into  (122)  leaves 


n 

X]  [amffm(a:.)  +  •  •  •+  aiffi(®i)  +  aoffo(*i)  -  /(®i)]  Si(®i)  =  0 

t  =  0 


for  j  =  0, 1,2, .  ..,m 


(123) 


Multiplying  through  on  the  left-hand  side 


n 

m9m{xi)9i{xi)  +  •  •  •  +  aigi{xi)gj{xi) 

«=o 


+ao5o(j:i)Pi(a:i)  -/(*»•  )yi(®»')]  =  0 


for;  =  0,1,2,  ...,m 


Breaking  up  the  left-hand  side  into  separate  sums  results  in: 


n 

'^^m9m{Xi)gj{Xi)  + 
«=0 

-I- 


n 

•••  +  X3ai</i(x,-)ff;(x,) 

1=0 


n 


n 


Y,aogo{xi)gj{xi)  -  =  0 

1=0  1=0 


for;  =  0,1,2, 


Next  factor  the  o’s  out  of  each  term: 


n 

®m  Xy  9m{xi)gj{Xi) 
1=0 


n 

+  ••■  + 01  X^ifi  (*•)»;■(*«•) 

{=0 

n  ti 

+  aoY^go{xi)gj{xi)  =  53/(x,)if;(x,) 

1=0  1=0 


(124) 


for;  =  0,l,2,---,m 
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Eq  (124)  provides  a  useful  result  because  each  of  the  summations  can  be  evaluated  since  the  gi's 
are  known  functions  and  the  x,-’s  are  tabulated  values.  Eq  (124)  is  actually  a  system  of  m+1  linear 
equations  with  the  Oj’s  the  only  unknowns.  These  linear  equations  can  then  be  solved  by  using  any 
of  a  number  of  standard  methods  for  the  solution  of  a  system  of  simultaneous  linear  equations. 

To  simplify  the  process  a  bit  and  make  it  more  suitable  for  computer  solution  it  should  be 
noted  that  the  coefficient  of  at  in  the  jth  equation  is 


okj  =  Yl3kixi)gj{xi) 


1=0 


for  iS:  =  0,1,2,..., m 
i  =  0,l,2 . m 


(125) 


Since  j  and  k  in  Eq  (125)  can  be  interchanged,  atj  =  ajk,  which  reduces  the  work  required  to 
compute  these  coefficients  by  half. 

The  method  presented  provides  a  workable  scheme  for  finding  the  least-squares  fit.  Using 
Eq  (124)  the  m  +  1  equations  are  solved  for  the  unknown  a’s,  which  exist  and  are  unique  if  the 
g(x)  functions  are  chosen  properly.  The  resulting  o’s  are  then  substituted  into  Eq  (112)  to  give  the 
desired  expression  for  pm(®)- 
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IV.  Results 


To  determine  the  surface  current  and  charge  distribution  of  the  antenna  by  using  the  method 
of  moments,  an  efficient  method  is  needed  to  numerically  evaluate  the  Green’s  functions.  The 
Green’s  functions  as  developed  by  Nazar  are  complex  valued  integral  functions.  In  his  thesis  he 
examined  the  Green’s  functions  expressions,  developed  techniques  to  numerically  integrate  them, 
and  coded  representative  examples  to  show  their  behavior  as  a  function  of  the  radial  separation 
between  the  source  and  observer.  This  thesis  effort  concentrated  on  coding  the  remaining  Green’s 
functions  needed  and  developing  approximations  that  could  be  used  in  an  efficient  moment  method 
analysis  of  the  antenna.  All  the  needed  Green’s  functions  were  calculated  numerically  and  then 
various  curve  fitting  schemes  were  tried  until  suitable  approximations  were  found  for  each.  The 
curve  fitting  was  done  using  a  least  squares  method  to  test  the  fit  of  the  various  approximations. 
All  the  custom  written  code  developed  in  this  thesis  was  written  in  FORTRAN. 

4.1  Numerical  Evaluaiion  of  the  Greenes  Funciions 

Before  approximations  could  be  made  for  the  Green’s  functions,  the  remaining  Green’s  func¬ 
tions  needed  to  be  coded  and  numerical  results  obtained.  The  examples  coded  by  Nazar  represented 
all  the  unique  forms  of  the  Green’s  functions  so  his  results  were  extended  by  coding  the  remaining 
Green’s  function  integrands  and  using  the  integration  routines  he  developed  [3].  A  complete  list  of 
the  Green’s  functions  needed  for  this  model  is  listed  in  Table  3.  The  redundant  forms  are  indicated 
by  their  equivalent  forms.  The  Green’s  functions  used  to  calculate  the  contribution  to  the  E  field 
from  the  equivalent  magnetic  sources  have  an  angular  as  well  as  a  radial  dependence.  Each  of  these 
functions  is  composed  of  two  or  three  integral  functions  which  each  have  only  a  radial  dependence. 
Table  4  shows  which  integral  equations  are  the  components  of  these  Green’s  functions.  All  the 
approximations  in  this  thesis  have  only  a  radial  dependence. 
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Table  3.  Green’s  Functions  for  Various  Source  Contributions 


HED  on  Interface  2b 

HED  on  Interface  3b 

G^’S'aCR) 

g^=S'3(R)=g^1,mR) 

G^33(R) 

G^“S(R)=G^’‘3='2(R) 

g^"2"3(R)=g^’^MR) 

HMD  on  Interface  lb 

G^aKR-O 

G|^j^(R.C) 

1 

g|^3’1(R.C) 

g|^3'i(R,C) 

G|^3"i(R,C) 

Table  4.  Green’s  Functions  Integral  Components 


HMD  on  Interface  lb 

G|*2l(rb2i(R).i*b2i(R)) 

G^^2l(I£2l(R).I*b2l(R).I‘b2l(R)) 

G^^2'la&2l(R),I'b2l(R).IUl(R)) 

G|f,{I^2l(R).I*b2l(R)) 

G^|1(IL3l(R).I*b3l(R)) 

G^l(Fb3i(R).ii3i(R).iL3i(R)) 

G||=i(iL3i(R).i|3i(R).4i(R)) 

G^^3",(I[3l(R).I*b3l(R)) 
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4.1.1  Behavior  as  R-*  0.  To  determine  the  currents  using  the  method  of  moments,  it  is 
necessary  to  perform  a  surface  integration  over  the  source  and  observer  areas  [3].  When  the  source 
and  observer  are  on  separate  planes,  the  Green’s  functions  integrands  decay  exponentially  and  the 
Green’s  functions  can  be  solved  and  evaluated  for  a  discrete  series  of  R  values  between  =  0  and 
the  maximum  radial  separation  expected  in  the  antenna  design.  Once  the  Green  s  functions  have 
been  evaluated  for  a  discrete  number  of  R  values,  then  a  suitable  approximation  csm  be  made  to 
fit  to  the  data.  This  approximation  would  be  much  more  computationally  efficient  than  trying  to 
evaluate  the  Green’s  function  each  time  a  value  is  needed.  Unfortunately  this  method  fails  when 
the  source  and  observer  are  on  the  same  plane. 

When  the  source  and  observer  are  on  the  same  plane  the  integrands  of  the  Green’s  functions 

do  not  decay.  Nazar  developed  asymptotic  approximations  for  all  the  Green’s  functions  where 

« 

the  source  and  observer  are  coplanar.  For  example,  the  asymptotic  approximation  for  G^xxCr)  as 
R-*0  becomes 


fin 

2n 


1  M  nN  X  'J2b  C0sh(U2b(b2b  “  bib))  +  /it23U3b  sinh(u2b(b2b  -  bib))  J  , 

io  - ^(A) - 

Jo(x) 


1  +  fibiS 


1  +  fibl2  +  fibiZ  +  ^13 


/ 


^dx  +  i- - - 

R  R 1  +  (Xii2  +  fibiz  +  fibiz 


(126) 


The  first  and  second  integrals  in  Eq  (126)  are  now  over  finite  intervals  of  integration  and  the 
solutions  finite  for  any  values  of  R.  If  the  Bessel  function  in  the  second  integration  is  expanded  into 
its  series  expansion  and  integrated  term  by  term,  the  value  of  the  integral  goes  to  Ac  as  i?  -♦  0. 
The  surface  integration  of  the  last  term  in  Eq  (126)  needs  to  be  done  analytically  for  the  method  of 
moments  analysis.  This  asymptotic  approach  is  only  needed  for  small  R  for  the  Green’s  functions 
where  the  source  and  observer  are  on  the  same  plane.  The  asymptotic  approximations  are  given 
along  with  the  general  form  of  the  Green’s  functions  in  Appendix  B. 

Depending  upon  the  Green’s  function,  the  contribution  from  the  last  term,  which  is  integrated 
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analytically,  will  be  entirely  real  or  imaginary.  Therefore  at  least  one  component  of  the  Green’s 
functions  where  the  source  and  observer  are  on  the  same  plane  can  be  approximated  for  all  values 
of  R  of  interest. 

4.J.S  Evaluation  Intervals  and  Antenna  Parameters.  All  the  unique  Green’s  functions  listed 
in  Table  3  were  numerically  evaluated  using  the  techniques  developed  by  Nazar.  These  functions 
were  the  ones  necessary  to  describe  the  tangential  electric  fields  on  both  antenna  patches.  All  the 
Green’s  functions  where  the  source  and  observer  are  on  the  same  plane  were  solved  for  50  values  of 
R  between  0.001  and  0.1  meters.  The  asymptotic  forms  of  these  functions  were  solved  for  25  values 
of  R  between  0  and  0.003  meters.  The  functions  where  the  source  and  observer  were  on  different 
planes  were  solved  for  50  values  of  R  between  0  and  0.1  meters.  A  quadratic  distribution  [17,  3]  of 
points  was  used  so  that  a  higher  concentration  of  data  points  would  be  taken  near  Rmin  where  the 
function  changes  more  rapidly.  This  quadratic  distribution  can  be  described  by 

Ri  =  O'  -  1)  +  Rmin  (127) 

where  N  is  the  total  number  of  points  to  evaluate.  The  functions  were  evaluated  for  the  following 
antenna  parameters,  which  were  the  same  as  used  by  Terry  in  his  work  [5]. 

Frequency  =  3.7Ghz 
6ij  =  0.00158m 
b2i  =  0.00316m 
cit  =  2.2eo 
£26  =  2.2£o 
£36  =  Eo 

=  A‘26  =  /J36  =  tio 


The  data  from  the  numerical  evaluation  of  the  Green’s  functions  was  then  approximated  using  the 
functions  described  in  the  next  section. 
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Green’s  Function  Approximations 


The  Green’s  functions  are  complex  valued  integral  functions.  The  purpose  of  this  thesis  was 
to  determine  approximations  for  the  real  and  imaginary  components  which  could  be  used  in  a 
method  of  moments  analysis.  Most  functions  can  not  be  evaluated  exactly  even  though  we  may 
handle  them  as  if  they  were  known  exactly.  All  functions  can  be  approximated  by  a  polynomial 
or  a  rational  function,  which  is  a  polynomial  divided  by  another  polynomial  [20].  At  the  outset  of 
this  effort  it  was  hoped  that  polynomial  approximations  could  be  used  for  all  the  components  of 
the  Green’s  functions,  but  this  was  not  the  case. 

All  of  the  Green’s  functions  analyzed  had  the  same  general  form.  They  are  complex  valued 
functions  involving  the  integration  of  a  Hankel  function  over  an  infinite  interval.  Nazar  used  a 
Hankel  transform  to  transform  all  the  Green’s  functions  to  integrals  of  a  Bessel  function  over  a 
semi-infinite  interval.  The  only  real  difference  in  the  functions  is  caused  by  the  multiplicative 
factor  in  front  of  the  integrand.  Some  were  real  and  some  were  imaginary.  The  Green’s  functions 
accounting  for  the  contribution  to  the  E  field  from  electric  and  magnetic  surface  currents  had  a 
real  factor  out  front,  while  those  accounting  for  the  contribution  to  the  E  field  from  electric  surface 
charge  had  an  imaginary  factor.  This  can  be  seen  from  derivation  of  the  electric  scalar  potential 
from  magnetic  vector  potential  by  using  the  Lorentz  gauge  condition  (see  Eq  (77)).  As  an  example 
G^’22(R)  and  Gq22(R)  are  shown 


nbxx 

'^AZZ 


(AR)A 


U2bCosh(u2b(b2b  -  bib))  -f  A^^23U3bsinh(u2b(b2b  -  bib)) 
D^(A) 


dA  (128) 
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rk  (m-  1  f 

'■■4.ju£2j.DS(i)DS,(A) 


e623U2b“3bCOsh2(u2b(b2b  “  bib)) 

/i423U^b“3b  +  iHbliVilb  +  A‘613e623U§b) 
♦uibtanh(bibUib) 
*sinh®(u2b(b2b  -  bib)) 

«2b  +  /'423e423U3b  +  i^bl3  +  /iH2C423) 

*uibU3btanh(bibUib) 


VdA  (129) 


♦ii§^sinh(2u2b(b2b  -  bib)) 

+/i4i3cmUibU^btanh(bibUib) 

V  / 


Notice  G^22(^)  bas  a  1/j  dependence  that  Ga22(^)  "^be  net  effect  of  the  1/j  factor  is  to 

interchange  the  real  and  imaginary  components  and  change  the  sign  of  the  real  component.  With 
these  differences  addressed  we  can  now  move  to  the  approximations. 


4.S.J  Single  Polynomial  Approximation.  At  the  outset  of  this  thesis  it  was  hoped  to  find  a 
polynomial  approximation  for  each  component  of  the  Green’s  functions  that  would  be  valid  over  the 
entire  range  of  radial  separations  of  interest.  This  was  only  possible  with  half  of  the  components. 
The  graph  of  the  imaginary  part  of  G^’22(R)  shown  in  Figure  10  illustrates  the  behavior  where  a 
single  polynomial  approximation  is  valid  over  the  entire  range  of  radial  separations.  The  behavior 
is  that  of  a  slowly  decaying  oscillatory  function  which  can  be  approximated  by  a  single,  sixth  order 
polynomial  of  the  form 


f(R)  —  ao  +  ai  R  +  a2R^  +  asR^  +  a4R‘*  +  asR®  +  asR®  (130) 


The  sixth  order  was  used  to  ensure  a  good  fit  for  all  the  components  of  the  Green’s  functions 
examined. 


The  approximations  were  chosen  basically  by  a  trial  and  error  method.  Several  trial  form 
polynomials  were  chosen  and  a  least  squares  curve  fitting  method  was  used  to  determine  the  co- 


Figure  10.  Plot  of  Imaginary  Part. 

efficients.  The  trial  polynomials  were  then  evaluated  at  the  same  values  of  radial  separation  as 
the  original  data.  Then  the  original  and  approximated  data  were  plotted  on  the  same  graph  and 
analyzed  to  see  how  well  they  agreed.  Several  iterations  were  enough  to  determine  the  order  of  the 
approximations  necessary  to  match  the  behavior  exhibited  in  Figure  10. 

4.2.S  Piecewise  Approximaiion.  The  behavior  of  the  other  half  of  the  components  to  be 
approximated  is  illustrated  by  the  graph  of  the  real  portion  of  0^22(11')  shown  in  Figure  11.  This 
behavior  is  that  of  a  sharply  decaying  oscillatory  function.  At  full  scale  the  oscillatory  nature 
is  not  evident  (see  Figure  11a),  but  if  the  scale  is  expanded  as  in  Figure  11b  it  can  be  seen. 
Numerous  attempts  were  made  to  approximate  curves  with  this  behavior  by  using  a  single  function. 
Polynomials  up  to  and  including  tenth  order  were  tried.  Good  agreement  was  attained  over  the 
sharply  declining  portion  but  the  oscillatory  portion  had  too  much  variation. 

The  interval  of  interest  was  broken  into  two  parts.  The  first  was  0.001-0.010  meters.  On 
this  interval  the  function  changed  the  most  rapidly  and  an  eighth  order  polynomial,  as  shown  in 
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Eq  (131),  was  used  to  approximate  all  the  Green’s  functions  over  this  interval.  As  before  the  least 
squares  curve  fitting  method  was  used  to  test  the  different  curves  and  determine  the  polynomial 
coefficients.  The  generality  of  using  a  single  function  for  the  approximation  was  lost  but  it  was 
possible  to  use  the  same  interval  breakpoint  for  all  the  Green’s  functions  approximated. 

g(x)  =  bo  +  biR+  baR^  +  baR®  +  b4R^  +  bjR®  +  boR®  +  brR^  +  bgR®  (131) 

Over  the  interval  from  0.01  to  0.1  meters  the  approximating  function  needed  to  match  the  oscillatory 
nature  as  well  as  the  remaining  portion  of  the  sharply  decaying  Green’s  functions.  For  an  arbitrary 
breakpoint,  a  simple  polynomial  would  not  match  the  curve.  Including  several  terms  with  negative 
powers  of  R  was  tried  in  order  to  match  the  sharply  decaying  portion.  These  curves  were  found 
to  be  well  approximated  by  an  eighth  order  polynomial  plus  two  terms  with  negative  powers  of  R. 
This  is  actually  equivalent  to  a  rational  function  of  two  polynomials.  The  numerator  would  be  a 
tenth  order  polynomial  and  the  denominator  would  be  a  quadratic  polynomial  with  the  constant 
and  first  order  term  equal  to  zero.  Rational  functions  of  polynomials  generally  give  slightly  more 
efficient  approximations  but  are  harder  to  obtain  [20].  The  coefficients  for  this  approximation  were 
also  obtained  by  using  the  least  squares  method  with  two  of  the  terms  having  negative  powers. 

h(x)  —  c_2R  ^  +  c_2R  ^  +  Co  +  ciR  +  C2R^  +  C3R®  +  C4R**  +  C5R®  +  c$R®  +  ctR^  +  cgR®  (132) 

A  listing  of  all  the  Green’s  functions/components  with  the  form  of  the  applicable  approximation 
for  each  range  is  given  in  Appendix  C. 

4.S  Sample  Results  of  the  Approximations 

To  show  how  well  the  approximations  matched  the  Green’s  functions,  sample  plots  of  the 
original  numerical  data  with  the  approximations  superimposed  are  displayed  in  the  next  several 


sections. 


4.3.1  Results  for  ^qZzC^)-  ®Axx(^)  Gq22(f^)  chosen  to 

illustrate  the  behavior  of  the  Green’s  functions  when  the  source  and  observer  are  on  the  same 
plane.  The  Green’s  functions  are  plotted  from  0.001  -  0.100  meters.  The  asymptotic  expressions 
were  not  plotted  along  with  the  full  expressions.  As  stated  in  Section  4.1.1,  one  component  of  the 
asymptotic  form  is  defined  at  iZ  =  0  while  the  other  goes  to  infinity  because  of  the  1/R  dependence 
in  the  third  term  of  the  asymptotic  form  illustrated  in  Eq  (126). 

Figure  12  is  a  plot  of  the  real  portion  of  GaxxCR))  the  original  and  approximated 
data.  Figure  12a  is  the  full  scale  data  and  Figure  12b  is  the  expanded  data  showing  the  decaying 
oscillatory  nature  of  the  Green’s  functions.  The  eighth  order  and  rational  polynomials  are  excellent 
approximations  for  the  real  portion  of  the  function.  The  imaginary  part  of  Ga^xx(^)  shown  in 
Figure  13.  The  much  more  slowly  decaying  oscillatory  nature  is  well  approximated  by  a  sixth  order 
polynomial.  The  asymptotic  form  of  the  imaginary  component  of  G^xx(R)  finite  at  iZ  =  0.  In 
fact  the  value  is  approximately  constant  over  the  range  from  0-0.001  meters. 

The  real  portion  of  Gq22(Il')  plotted  in  Figure  14  shows  the  function  is  well  approximated 
by  a  sixth  order  polynomial.  The  value  of  the  asymptotic  form  is  almost  constant  over  the  range 
0-0.001  meters.  This  behavior  can  be  seen  in  Figure  14  for  values  of  radial  separation  close  to  zero. 
The  plot  is  of  the  same  shape  as  the  imaginary  portion  of  Gaxx(^)-  switch  in  the  behavior 
of  the  real  and  imaginary  components,  as  compared  to  G^xxC^)*  Vi  factor  in  the 

Green’s  function. 

The  imaginary  portion  of  Gq22(Il')  is  plotted  in  Figure  15  and  is  of  similar  shape  to  the  real 
part  of  Gaxx(^)  ®x®®Pt  for  ffi®  of  the  values.  Here  again  the  eighth  order  polynomial  and  the 
rational  polynomial  provided  excellent  approximations  for  the  original  function. 
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Figure  12.  Plot  of  Ga22(R')  Rea!  Part, (a)  Full  Scale,  (b)  Expanded  Scale. 
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Figure  15.  Plot  of  Imaginary  Part,(a)Full  Scale,  (b)  Expanded  Scale. 


4.5.2  Results  for  GqjaW-  ®Axx(^)  ®q23(^)  chosen  to 

illustrate  the  behavior  of  the  Green’s  functions  when  the  source  and  observer  ate  on  different 
planes.  These  functions  are  plotted  from  0  -  0.100  meters  and  no  asymptotic  forms  were  needed. 

Figure  16  is  a  plot  of  the  real  portion  of  original  and  approximated 

data.  Figure  16a  is  the  full  scale  data  and  Figure  16b  is  the  expanded  data  showing  the  decaying 
oscillatory  nature  of  the  Green’s  functions.  The  eighth  order  and  rational  polynomials  are  again 
excellent  approximations  for  the  real  portion  function. 

The  imaginary  part  of  G^xx(^)  shown  in  Figure  17  is  of  similar  shape  as  the  imaginary  part 
of  G^xx(^)‘  much  more  slowly  decaying  oscillatory  nature  is  again  well  approximated  by  a 
sixth  order  polynomial. 

The  real  portion  of  Gq23(R)  is  plotted  in  Figure  18.  This  function  is  well  approximated  by 
the  sixth  order  polynomial.  The  plot  is  of  the  same  shape  as  the  real  portion  of  Gq22(R)'  Again 
the  switch  is  due  to  the  1/j  factor  in  the  Green’s  function. 

The  imaginary  portion  of  Gq23(R)  is  plotted  in  Figure  19  and  is  of  similar  shape  to  the 
imaginary  part  of  Gq22(R)-  Here  again  the  eighth  order  polynomial  and  the  rational  polynomial 
provided  excellent  approximations  for  the  original. 

4.5.3  How  Well  the  Approximations  Matched.  The  initial  tests  as  to  how  well  the  approxi¬ 
mations  matched  the  Green’s  functions  were  done  by  “eye”.  The  approximations  for  all  the  Green’s 
functions  basically  overlaid  the  graph  of  the  functions  they  were  to  approximate.  As  a  further  test 
of  the  “goodness”  of  fit,  each  function  and  its  approximation  was  integrated  over  the  entire  range  of 
values  of  R  of  interest  using  a  simple  trapezoidal  integration  routine.  The  results  of  these  integra¬ 
tions  for  the  functions  displayed  in  Sections  (4.3.1)  and  (4.3.2)  are  given  in  Table  5.  The  agreement 
between  the  actual  values  and  the  approximations  was  excellent.  Comparing  the  actual  and  ap¬ 
proximated  results  show  agreement  out  to  the  third  or  fourth  decimal  place,  indicating  agreement 
to  within  a  fraction  of  a  percentage. 
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Figure  16.  Plot  of  G^23(^)  Combined  Data  (a)Full  Scale,  (b)  Expanded  Scale. 
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0.1  so 


Figure  17.  Plot  of  G^2’3(R)  Imaginary  Part,  Combined  Data  . 


Figure  18.  Plot  of  Gq23(R)  Real  Part,  Combined  Data. 
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Figure  19.  Plot  of  Gq23(R)  Imaginary  Part,  Combined  Data  (a)Full  Scale,  (b)  Expanded  Scale. 
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Table  5.  Integration  of  the  Green’s  Functions  and  their  Approximations 


Component 

Real 

Actual 

Approx 

Actual 

Approx 

8.2416e-8 

8.2426e-8 

-2.3774e-9 

-2.3774e-9 

G522(R) 

3.5832e-8 

3.5817e-8 

-1.5358e-l 

-1.5359e-l 

G^=‘2l.(R) 

1.1871e-7 

1.1871e-7 

-4.8477e-9 

-4.8477e-9 

7.5807e-3 

7.5776e-3 

-2.5025e-l 

-2.5025e-l 

Data  Supplied  to  Capt  Irvin 

The  coefficients  for  the  approximations  of  the  full  expression  Green’s  functions  were  given  to 
Capt  Irvin  for  use  in  his  moment  method  analysis  of  the  antenna.  A  complete  set  of  the  coefficients 
is  given  in  Appendix  C.  The  data  for  the  asymptotic  forms  was  given  to  him  in  tabular  form  and 
he  interpolated  to  find  any  desired  value.  The  handling  of  the  asymptotic  data  was  dictated  by 
time  constraints  and  by  some  uncertainty  at  the  time  about  evaluating  the  asymptotic  forms  of  the 
Green’s  functions  at  i?  =  0.  As  stated  in  Section  (4.1.1),  one  component  of  the  approximation  could 
be  incorporated  into  the  full  expression  because  it  if.  finite  at  i?  =  0.  The  term  in  the  asymptotic 
approximation  with  the  l/R  dependence  dominates  the  behavior  of  the  component  which  is  not 
finite  at  iZ  =  0.  Perhaps  the  asymptotic  data  for  small  values  of  ii  0  could  be  incorporated  into 
the  full  expression  approximation.  This  would  only  leave  the  analytical  surface  integration  of  the 
last  term  at  =  0. 
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V.  Conclusions 


Numerical  approximations  for  a  theoretical  model  of  an  aperture  fed  stacked-patch  microstrip 
antenna  have  been  developed.  The  behavior  of  the  Green’s  functions  were  examined  and  three 
different  approximation  forms  were  found  applicable  to  each.  The  approximations  were  then  used 
in  a  separate  thesis  effort  to  determine  the  patch  surface  currents  and  charge  distributions.  The  rest 
of  this  chapter  discusses  the  answers  to  the  research  questions  posed  in  the  introduction,  suggestions 
for  further  work,  and  some  general  observations  made  during  the  study. 

5.1  Answers  to  Research  Questions 

5.1.1  What  is  the  form  of  the  functions  necessary  to  accurately  model  the  Green’s  function? 
Each  Green’s  function  can  be  modeled  by  three  separate  approximations.  One  component  of  each 
Green’s  function  can  be  modeled  with  a  sixth  order  polynomial  over  the  entire  range  of  radial 
separations.  The  other  component  required  a  piece-wise  approximation  over  two  separate  intervals 
which  were  the  same  for  each  function.  Over  the  interval  from  Rmin  out  to  0.01  meters,  the  rapidly 
changing  functions  can  be  modeled  by  an  eighth  order  polynomial.  From  0.01  -  0.10  meters  an 
eighth  order  polynomial  plus  two  terms  with  negative  powers  were  used  to  model  the  functions. 
This  expression  is  actually  equivalent  to  a  ratio  of  polynomials  with  a  tenth  order  polynomial  in 
the  numerator  and  a  quadratic  polynomial  in  the  denominator. 

5.1.2  How  well  do  these  functions  approximate  the  original  Green’s  functions?  The  approx¬ 
imations  were  found  to  be  excellent  approximations  for  the  original  Green’s  functions.  When 
initially  developed  the  “goodness”  of  fit  was  determined  by  “eye”  and  found  that  the  approxima¬ 
tions  overlaid  the  original  functions  over  the  entire  range  of  radial  separations  values  of  interest. 
As  a  further  test  of  the  agreement  the  original  functions  and  the  approximations  were  integrated 
using  a  trapezoidal  integration  routine  and  the  results  were  found  to  agree  to  within  at  least  two 
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and  sometimes  and  possibly  four  decimal  places.  This  indicates  an  error  on  the  order  of  a  fraction 
of  a  percent. 

5.1.S  Can  the  approximations  be  simplified  while  still  retaining  the  accuracy  required  for  the 
model?  This  question  was  not  answered  due  to  lack  of  time  for  an  iterative  process.  The  actual 
approximations  were  determined  late  in  the  thesis  process  and  needed  by  Captain  Irvin  to  complete 
his  thesis  effort.  Continued  work  on  this  question  is  needed  and  discussed  in  the  next  section. 

5.S  Recommended  Follow  on  Work 

As  mentioned  in  the  previous  section,  work  is  needed  on  determining  if  simpler  approximations 
would  be  acceptable  in  the  model.  Further  work  on  these  approximations  needs  to  be  done  in  concert 
with  the  moment  method  work  done  by  Capt  Irvin.  Using  the  approximations  developed  in  this 
thesis  and  the  surface  currents  determined  by  Capt  Irvin  as  a  baseline,  the  effect  of  simplifying 
the  approximations  can  be  studied.  As  the  radial  separation  between  the  source  and  observer 
increases,  the  values  of  the  Green’s  functions  decrease  by  several  orders  of  magnitude.  The  effect 
of  determining  a  lower  threshold  value  of  the  Green’s  functions,  below  which  they  may  be  assumed 
to  be  zero  should  be  investigated. 

Another  area  of  interest  would  be  to  take  the  integration  routines  written  by  Nazar  and  try 
to  simplify  them.  The  IMSL  routines  used  in  his  code  are  extremely  elaborate.  Simpler  routines 
may  speed  up  the  numerical  evaluations  of  the  Green’s  functions  required  for  each  set  of  antenna 
parameters.  For  example  the  execution  time  of  the  modified  Nazar  code  used  in  this  thesis  was  one 
hour  on  an  ELXSI  System  6400  Computer,  running  under  its  native  EMBOS  operating  system. 

The  code  developed  in  this  thesis  needs  to  be  combined  into  one  routine  to  determine  the 
Green’s  functions  approximations.  Separate  routines  were  written  for  each  approximation  procedure 
but  the  data  input  was  modified  to  include  only  the  applicable  portions.  It  should  be  a  simple 
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matter  to  incorporate  the  routines  into  one  program  with  conditionals  to  determine  the  appropriate 
approximations  to  use. 


61 


Appendix  A.  Vector  Potential  Parameters 


A.l  HMD  on  Interface  lb. 


N^xi(A) 

N^xi(A) 

N^xi(A) 


=  eijA 


NLW  = 


— eijA  [ei23U3bCosh(b2bU2b)  +  U2bsinh(b2bU2b)]  sech(bibUib) 
eijA  [u2bCosh(b2bU2b)  +  e423U3bsinh(b2bU2b)]sech(bibUib) 
ei6Au2Bexp(b2bU3b)sech(bibUib) 

[uib  +  e»i3U3btanh(bibUib)]  U2bCOsh(u2b(b2b  -  bib)) 

+  [£423UlbU3b  +  £412U^b*'a”^(blbUlb)]  sinh(U2b(b2b  “  bib)) 
[-M413£413  +  COsh^(U2b(b2b  “  bib))]  2u^bCsch(2bibUib) 

+  [-P412£412uib  -  Wl3£413«3b  +  /^l23e423U§b] 
*2csch(2bibUib)sinh^(u2b(b2b  -  bib)) 

+  [/^423  -  /i413e412  “  /i412£413  +  £423] 

*U2bU3bCsch(2bibUib)sinh(2u2b(b2b  -  bib)) 

[1  -  /i423£423]£i4UibU2bCsch(bibUib)cosh(bibU2b) 

[-^'412  +  /<4i3£423]£i4U2bSech(bibUib)sinb(bibU2b) 

(A»413£14  -  /i423£24)U3bCOsh(b2bU2b) 

+(|i412£l4  -  £24)U2bSinh(b2bU2b) 
*U2bsech(bibUib)cosh(u2b(b2b  -  bib)) 

(;i413£413  -  A‘423£423)U3bCOsh(b2bU2b) 

+(/i412£413  -  £423)u2bSinh(b2bU2b) 
*£24U3bsech(bibUib)sinh(u2b(b2b  -  bib)) 


+ 
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Nax2(-^)  =  [u2bCosh(u2b(b2b  -  bib))  +  ;^i23U3bsinh(u2b(b2b  -  bib))]  csch(bibUib) 

Nbx2(^)  =  -/il6-^[W23U3bCOsh(b2bU2b)  +  U2bSinh(b2bU2b)] 

Ncx2(-^)  =  A‘U^[U2bCOsh(b2bU2b)  +  /'M3«3bSinh(b2bU2b)] 

Ndx2(^)  =  /^li^exp(b2bU3b)U2b 
» 

[1  -  /l613eM3]u2bS®ch(bibUib) 

+  [(1  -  /iM2e612)uib  +  i-mi3^bl3  +  /i623E623)uib] 

Nb22('^)  =  '  ♦sech(bibUib)sinh^(u2b(b2b  -  bib)) 

+[/^t23  -  tlh\3^b\2  —  lib\2^b\3  +  £623] 

♦i^2i^sech(bibUib)sinh(2u2b(b2b  -  bib)) 
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N^z2(A)  = 


N^z2(A)  = 


Nd22(A)  =  exp(b2bU3b)u2bA^  < 


[-£H2  +  A‘i23en3]/ii»u^bSin5i(bibU2b) 

[1  -  A‘623e623]A‘i6UibU2bCOsh(bibU2b)tanh(bibUib) 

(Miie6i3  -  A‘2l£m)u3bCosh(b2bU2b) 

+(-^2i.  +  /iiie»i2)u2bsinh(b2bU2b) 

♦U2bCOsh(U2b(b2b  -  bib)) 

ifibiz^bia  —  M623e»23U3bCOsh(b2bU2b) 

+(-/^623  +  Wl3eH2)U2bSinh(b2bU2b) 
*/i26U3bsinh(u2b(b2b  -  bib)) 

(e»12  -  A‘623eM3]/lUU2b“sh(bibU2b) 

[-1  +  /i623e»23]A'uUibU2bsmh(bibU2b)tanh(bibUib) 

(-/iiieM3  +  M2(.ek23)u3bsinh(b2bU2b) 

+{l^2b  -  /iUeM2)U2bCOSh(b2bU2b) 

*U2bCOsh(U2b(b2b  -  bib)) 

(“A'Jiacua  +  A'623em)u3bsinh(b2bU2b) 

+(/i»23  -  /i613eM2)U2bCOSh(b2bU2b) 
*/i2lU3bsinh(u2b(b2b  -  bib)) 

[a‘36  -  A»iteii3]u2bCOsh(u2b(b2b  -  bib)) 
(/i26  -AilieH2)U3b 

•i-{li3blibi2  -  /iue»23)uibtanh(bibUib) 

♦sinh(u2b(b2b  -  bib)) 


A.S  HED  on  Interface  3b. 


N^aCA)  =  /iuAu2bCSch(bibUib) 


64 


Nbx3(^)  =  /i26A[uibCOsh(bibU2b)coth(bibUib)  -  ^ii,i2U2bsinli(bibU2b)] 
NcxsC^)  =  /i2iA[-uibsinh(bibU2b)coth(bibUib)  +  /iH2U2bCosh(bibU2b)] 


Ndx3(^)  =  /<2lAexp(b2bU3b) 


/iil2U2bCOsh(U2b(b2b  -  bib)) 
+UibCoth(bibUib)sinh(u2b(b2b  -  bib)) 


[1  -  /i6i3CM3]u2bCosh(u2b(b2b  “  bib))sech(bibUib) 


^  (efcl2  “  /i^23CM3)WlbCSCh(bibUib) 

+ 

+(-/^ki2e6i3  +  em)u3bsech(bibUib) 


sinh(u2b(b2b  -  bib)) 


NL(A)  = 


Nc"z3(A)  = 


(/iueil3  -  /i26e623)/^3iCOSh(b2bU2b) 

U2b 

+(-/i2i  +  /iueM2)/i26sinh(b2bU2b) 

(-en2  +  /i423en3)u2bsinh(bibU2b) 

+ 

+(1  -  /i423e»23)uibCosh(bibU2b)tanh(bibUib) 
♦/illU2bCOsh(U2b(b2b  -  bib)) 

(1  -  /i623e423)uibCOSh(bibU2b) 

+ 

+(-C412  +  /X423e413)U2bCOth(bibUib)sinh(bibU2b) 
*/i24Uibsinh(u2b(b2b  -  bib)) 

{-IJtihCbia  +  /i24e423)#i34sinh(b2bU2b) 

U2b 

+(+/i26  -  A<16e612)M26COSh(b2bU2b) 

(ejl2  -  /il23eM3)U2bCOsh(bibU2b) 

+ 

+(-l  +  /i423e»23)uibsmh(bibU2b)tanh(bibUib) 
*A»uU2bCosh(u2b(b2b  -  bib)) 

(—1  +  /i»23e423)uibsinh(bibU2b) 

+ 

+(£412  -  ^423£413)U2bCOth(bibUib)cOsh(bibU2b) 
*M24Uibsinh(u2b(b2b  -  bib)) 
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Ndz3(^)  =  exp(b2bU3b)A*  < 


[-ufb+uy/i3i 

+  [/isfcuJb  -  /iucuauib]  cosh2(u2b(b2b  -  bib)) 
[-/«2i£»23Uib  +  sinh*(u2b(b2b  -  bib)) 

(fi3b£hii  -  /i2keii3)coth(bibUib) 

+(/i3»/i4i2  -  Mi»c»23)tanh(bibUib) 

,^uii^gjjjjj^2u2b(b2b  —  bib)) 
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Appendix  B.  Complete  Gieen^s  Functions 


Db(A) 


[/^il2U3b  +  UibCOth(bibUib)]  U2bC08h(U2b(b2b  -  bib)) 

+  [^‘il2U2b  +  /‘i23UlbU3bCOt(bibUib)]  8inh(U2b(b2b  “  bib)) 

[cn3U3b  +  Uibtanh(bibUib)]  U2bC08h(u2b(b2b  -  bib)) 

+  [c4i2U2b  +  «»23UibU3btan(bibUib)]  sinh(u2b(b2b  -  bib)) 


Uib 


BJ  Green’s  Functions  for  Electric  Fields  From  Electric  Sources 
B.1.1  BED  at  Interface  Sb. 


gS(R) 


£l* 

47r 


= 


y^B(^)(l,p)\^jbCOsh(U2b(b2b  -  bib))  +  /^423U3bSinh(U2b(b2b  -  ^Ib)) 

( 

e423U2b«3bCOsh2(U2b(b2b  -  bib)) 

PmU2bU3b  +  (/ilnUjb  +  /JM3e423u|b) 
*uibtanh(bibuib) 

1  r  H^^^(AR)A  I  ♦8inh2(u2b(b2b  -  bib)) 

45rjwe2J  X  Db(A)Db  (A)  | 


+ 


'*2b  +  W23et23U3b  +  +  #14128423) 

*UibU3btanh(bibUib) 
*i*f^sinh(2u2b(b2b  -  bib)) 
+Ai4i3emuibU2btanh(bibUib) 
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For  large  A,  these  integrands  decay  as  A 


G») 


=  Gb-(R) 


H^^^(AR)Au2b 


[U3b  +  /iM3Uibtanh(bibUib)] 
*U2bCOSh(U2b(b2b  -  bib)) 

/‘ll3eH2U2b  +  (^^23  “  A‘tl3£H2)U3b 
+PM2UibU3btanh(bibUib) 

♦sinh(u2b(b2b  -  bib)) 


MA 


For  large  A,  these  integrands  decay  as  exp[— A(b2b  —  bib)]. 


B. 1.1.1  Asympioiic  Forms  for  Small  R. 


Gb-(R)  = 


G^22(R)  = 


^  i  r  Ia22(R.  '^)dA  -  7T - 

2’’'  1  7o  ^  ^  +  /'HZ  +  A'>23  +  /^H3  , 


Jo  ^ 


+ 


1 


1  +  PJ23 


R.  1  +  t^bl2  +  P»23  +  f^hl3 


} 


1 


r^c 


1  +  Si23 


27rjwe2i  (70 
1  1  +  £623 


+ 


R.  1  +  £612  +  £623  +  £613 


I‘'^^(R’A)dA-r+eM2  +  £623  +  £6J3 

} 


yRAc 

Jo 


Jo(x) 

R 


dx 


where  Ia22(R)  A)  and  Iq22(R'i  A)  are  the  integrands  of  the  original  functions  with  Ho^^(AR)  replaced 
by  Jo(AR). 


B.l.S  HED  ai  Interface  Sb. 


G^KR) 

g51^3"3(R)  =  G^^g^aCR) 


^  f  H^=^>(AR)A 

D^(A) 


471  ic 


P612U2bCOsh(u2b(b2b  “  bib)) 
+UibCOt(bibUib)sinh(u2b(b2b  -  bib)) 


dA 
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G533(R) 


1  f  h^^\ar)a 

4jrj«e36  Jc  D^(A)DS,(A) 


U3b  +  A‘H3Uibtanh(bibUib) 

♦cosh*(u2b(b2b  -  bib)) 

(uJb  +  A<M2eM2U2b)  «3» 

+A»l23eH2Ulbu|bCOth(blb^lb) 

♦sinh^(u2b(b2b  -  bib)) 


/l623Uib  +  MM3eH2U2b 
eM2Coth(bibUib) 


«2l 


+ 


+/Hi2tanh(bibulb) 


UibUsb 


J 


*J!“.sinh(2u2b(b2b  -  bib)) 


ua 


For  large  A,  these  integrands  decay  as 


G^»)  =  G^y/3(R)  =  G^f2(R) 


rb  rm  -  ^  /  H[,^>(AR)Au2b  I 

-  4,rwe2»  ic  D.^A)Db  (A)  1 


[(u3b  +  /iH3Uibtanh(bibUib)] 
♦ei23U2bCOsh(U2b(b2b  -  bib)) 

(-l+W23et23)u?b  +  uib 
+/iH2et23UibU3btanh(bibUib) 

*sinh(u2b(b2b  -  bib)) 


Ua 


For  large  A,  these  integrands  decay  as  exp[— A(b2b  —  bib)]. 


B. 1.2.1  Asympioiic  Forms  for  Small  R. 


can) 


=  —It 

27r  \yo 


+ 


I^33(R,A)dA- 

1  +  Wl2 


1  +  MH2 


1  +  fihl2  +  l^b2S  +  f^hl3 


,RAc 

Jo 


Jo(x) 

R 


dx 


R  1  +  t^H2  +  t^b2S  +  I^bl3 . 
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gUR)  = 


— (r 

2jrjwe3k 


l533(R,A)dA- 


l  +  gH2  Jo(x) 

1  +  SH2  +  £123  +  £M3  Jo 


dx 


2. _ 1  +  gtl2 _ 1 

R  1  +  £612  +  £623  +  £613  / 


where  Ia33(Ri  A)  and  Iq33(R,  A)  are  the  integrands  of  the  original  functions  with  H^*^(AR)  replaced 
by  Jo(AR). 

B.2  Green's  Funclions  for  Electric  Fields  From  Magnetic  Sources 
B.S.J  HMD  at  Interface  lb. 


GEfi(R.O  = 

G|"2l(R.O  = 

G^7i(R,0  = 
G^^2^(R.C)  = 


-sin(20  [5lC2i(R)-ll2i(R) 

-  l-Ib2i(R)  -  cos2(C)rb2i(R)  +  cos(2C)rb2,(R)] 

-  [I^2i(R)  +  sin2(OI^2i(R)  +  cos(2C)It2i(R)] 

-sin(2C)  [-|l^2i(R)  +  I*b2i(R) 


where 


Ib2l(R)  = 


^  H^o^^(AR)A3sech(bibUib) 


De"(A)D^(A) 

[1  -  /I613e613]u^b 

+  [(1  -  /^612e6l2)U2b  +  (-/^613£613  +  /I623£623)u3b] 
♦Sinh2(u2b(b2b  -  bib)) 

'l‘(/^623  -  M613C6I2  "  /^612£613  +  £623) 

*i!^sinh(2u2b(b2b-bib)) 


}dA 
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f  H^i^^(AR)A^sech(bibUib) 

X"  RD^(A)D^(A) 

{1  -  /iiiseualuib 

+  [(1  —  /^»t2EH2)U2b  +  (“/^M3£H3  +  /i623£623)U3b] 

♦  '  +sinh2(u2b(b2b  -  bjb)) 

+(/423  -  /^413e612  -  /i612eM3  +  £423] 


*!i2i^sinh(2u2b(b2b  -  bib)) 


Ib2l(R) 


H^o®)(AR)Au2bsech(bibUib)  en3U3bCosh(u2b(b2b  -  bib)) 

nb  f  \\ 

+£4i2U2bfinh(u2b(b2b  -  bib)) 


For  large  A,  these  integrands  decay  as  exp(“Abib). 


C3||1(R.C)  = 


-sin(2C)  ilL3i(H-)-I'b3i(R) 


GesKR.  0  =  -  [-iUi(R)  -  cos2(c)rb3i(R) + cos(2c)rb3i(R)] 
G^7i(R,0  =  -[lt3i(R)  +  sin2(0rb3i(R)  +  cos(2C)I‘b3i(R)] 
g|^3i(R.C)  =  -sin(2C)[-il[3i(R)  +  I|3i(R) 


Ib3i(R)  = 


h[,^)(AR)A3 

D."(A)Db(A) 


[1  -  /r4i3£4i3]u^bSech(bibUib)cosh(u2b(b2b  -  bib)) 
(£412  -  A‘423£4l3)UlbU2bCSch(bibUib) 

+ 

+(-^412£413  +  £423)U2bU3bSech(bibUib) 

♦sinh(u2b(b2b  -  bib)) 
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ItaiCR) 


4i(R) 


[1  -  /'M3eii3luLsech(bibUib)cosh(u2b(b2b  -  bib)) 
(£412  -  ^423e413)uibU2bCSch(bibUib) 
+(-/<612£413  +  £423)U2bU3bSech(bibUib) 

^  *sinh(u2b(b2b  -  bib)) 

°  [g4i3U3bsech(bibUib)]dA 


L 


RD^(A)D^(A) 


> 


For  large  A,  these  integrands  decay  as  exp(-Abib). 
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Appendix  C.  Summary  of  Results 


Green’s  Function  Approximation  Forms 


g(R)  =  bo  +  biR  +  baR^  +  baR^  +  b^RHbsR^  +  beR^  +  bTR^  +  bgR® 

h(R)  —  c_iR  ^  +  c-aR  ^  +  Co  +  ciR+ caR^  +  caR^  +  C4R'^  +  C5R*  +  c$R®  +  C7R^  +  cgR® 

f(R)  —  sq  ajR  +  aaR^  +  aaR^  +  a4R^  +  asR®  +  asR® 


Table  6.  Coefficient  Order  for  the  Piecewise  Approximation 


Eight 

1  Order 

bo 

bi 

62 

63 

64 

bs 

be 

67 

bi 

Ratio 

C-2 

c-l 

Co 

Cl 

C2 

C3 

C4 

C5 

c$ 

C7 

cs 

Table  7.  Coefficient  Order  of  the  Single  Function  Approximation 


Sixth  Order 

Co 

ai 

02 

aa 

04  1 

as 

as 
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Coefficients  for  Each  Green’s  Function  Approximation 


Real  Component 


8th  Order  0.001-0.010 

2.82713S7E-04  -3.a8466E6E-01 
-4.2391776E+C9  3.9943867E+11 

Ratio  p. 010-0. 100 

1.1336337E-10  -7.7606745E-09 
1.2917902E-10  3.4S24106E-01 

1.7830531E!-03 


2.728S907E+02  -1.0962870E+05 
-2.0833405E+13  4.6074276E+14 

2.877B042E-07  -7.9761966E-'14 
-9.6808158E+00  1 .2043431E+02 


Inagin2u:y  Component 


6th  Order  0.001-0.100 

-8.2697879E-08  3.0039B34E-08  4.4481S33E-0B  3.7821366E-04 

2.3993749E-01  -7.4304B13E-01 

Real  Component 


8th  Order  0.001-0.010 

2.989B392E-04  -3.9713808E-01  2.7003667E+02 

-4.2349631E+09  3.9979300E+11  -2.0878692E+13 

Ratio  0.010-0.100 

2.6798747E-10  -7 . B940224E-09  1.8608318E-07 

-3.1416726E-11  3.4B344B6E-01  -8.3B78282E+00 

9.774BBB0E+02 


-1.0882422E'f0B 

4.621BB27E+14 

1.B038343E-14 

8.86B8B36E-)-01 


Imaginary  Component 
6th  Order  0.001-0,100 

-3.219376BE-07  1.1728108E-07  1.73B9691E-04  1.4780323E-03 

9.38413B4E-01  -2 . 90732S2E+00 

G^”(R) 


Real  Component 


8th  Order  0.000-0.010 

4.360B698E-OS  1.6428BdlE-03 

4,9368870E+08  -4,6983776E+10 
Ratio  0.010-0.100 

1 , 7029449E-10  -8 . 139924BE-09 
6.3114878E-11  3.4697986E-01 

1.4970486E+03 


-1.9408606C«-01 

2.4369181E-«-12 

2.66221S7E-07 

-9.2439316E+00 


1.0921B13E+04 

-B.31378B4E+13 

-4.1338B61E-14 

1.0942193E-f02 


2.7278B77E+07 

-2.8113896E-04 

-7.3386638E+02 

-2.3998044E-02 

2.7177121E+07 

-3.224288BE-04 

-4.628BBSBE-I-02 

-9.3812439E-02 

-3.0476872E+06 

-2.996327BE-04 

-6.3847649E+02 
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% 
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Imaginary  Component 
6th  Order  0.001-0.100 


-1.6312406E-07  4.7418351E-08 

4.6971681E-01  -1 .4547107E+00 

IL2i(R) 

Real  Component 

8th  Order  0.000-0.010 

8.8848203E-05 

7.1341936E-04 

-4.6862616E-02 

-2.3690010E+04  -5. 1497773E+03 
6.0581647E+16  -6.9634200E+18 
Ratio  0.010-0.100 

2.0098000E-t-09 

4.0033S36E'f20 

1.1664054E-(-ll 

-9.2918335E+21 

-2.4628222E+14 

-2.7567323E-01  4. 1969248E+01 

-1.7263071E+03 

6.7013484E-04 

2.5393666E+06 

-1.0559782E+00  -3.8054S23E+09 
-1.8720S67E+13 

Imaginary  Component 

6th  Order  0.000-0.100 

1.0860188E+11 

-1.3359528E+12 

7.9367587E+12 

7.3622798E+02  1 .6413327E+02 

-8.2057805E+09  2.9878758E+10 

Real  Component 

8th  Order  0.000-0.010 

-9.4660002E+05 

-S.S624401E+06 

6.6599743E+08 

-7.1631702E+04  -1 .6126649E+06 
1.1322934E+17  -2.2354787E+19 
Ratio  0.010-0.100 

1.6138682E+10 

1.6596165E+21 

-4.6090776E+12 

-4.5322535E+22 

2.0514026E+14 

-3.9145822E-01  6.0771449E+01 

-2.4640424E+03 

1.0734826E-03 

4.2372332E+06 

-1.6926680E+00  -6.6280178E+09 
-3.1213003E+13 

Imagineiry  Component 

6th  Order  0.000-0.100 

1.8569636E+11 

-2.2676061E+12 

1.3351982E+13 

1.4578382E+03  3.2605113E+02 

-1.6258111E+10  5.9202168E+10 

It2l(R) 

Real  Component 

8th  Order  0.000-0.010 

-1.8762160E+06 

-1.1014941E+07 

1.3194276E+09 

-1.1845706E+04  6.1244112E+03 

9.6364022E+15  -9.8801478E+17 

4. 872061 8E+08 
6.1860820E+19 

3.fl07918E+10 

-1.1172704E+21 

-4.6628645E+13 
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Ratio  0.010-0.100 

-1.4338001E-01  -8.8320753E+C0  4 . 8974870E+02  -1 .9777346E-04  -1 .6920208E+05 
3.2111913E-01  2.6480174E+08  -9 . 6649951E+09  1 .47297S5E+11  -1 .0371150E+12 
2.7871008E+12 

Imaginary  Component 

6th  Order  0.000-0.100 


3.6839408E+02  -1.0S65021E+02 
-1.2791397E+09  4.0846293E+09 

Real  Component 

8th  Order  0.000-0.010 

-2.1992704E+05 

-1.7804531E+06 

1.2335704E+08 

-3.5845677E+04  -1.8303401E+06 
6.8048203E+16  -7.7186016E+18 
Ratio  0.010-0,100 

3.3506986E-^09 

4.8748677E+20 

-4.0133292E+11 

-1.2131994E+22 

-1.6475295E+14 

-3,2220387E-01  -1 . 1195931E+01 

6.7963919E+02 

-2.7539833E-04 

-7.0817217E+04 

4.4877070E-01  1 .6953818E+08 

3.3469197E+12 

Imaginary  Component 

6th  Order  0.000-0.100 

-8.6162929E+09 

1.S436197E+11 

-1 . 1848456E+12 

7.2947448E+02  -2.0927389E+02 
-2.5345126E+09  8.0937999E+09 

Real  Component 

8th  Order  0.000-0.010 

-4.3564658E+0B 

-3.5273908E+06 

2.4440619E+08 

8.3730489E+05  1, 854891 lE+07 

3.0643030E+19  -3.2457653E+21 
Ratio  0.010-0.100 

-6.9596640E+11 

1.8334913E+23 

4.9585244E+14 

-4,2836938E+24 

-1.65090? 

1.0488087E+00  -9,7803738E+01 

3.8139458E-(-03 

-l,1916493E-03 

-4.520323iEvO 

1.8826109E+00  6.2316114E+09 

-1,7672322E+11 

2.1844777E+12 

-1.31£j428E+-'3 

3.1297587E+13 
Imaginary  Component 
6th  Order  0.000-0.100 

-1.2618585E+03  -1 . 1645510E+01  1 .2324177E+06  7.8552878E+06  -8. 1658149E+08 

9.7146319E+09  -3.4550352E+10 
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4 


IL3i(R) 


Real  Component 

8th  Order  0.000-0.010 


2.8114574E+06  6.9078766E+06 

-4.0848676E417  8.1118069E+19 

Ratio  0.010-0.100 

-5.9241731E+10 

-6.0339466E+21 

1.6899417E+13 

1.649E873E+23 

-7.73EE313E+14 

1.5538455E400  -1.2252199E402 

4.E582187E+03 

-1.4436E00E-03 

-6.4603264E+06 

2 . 2647814E400  9 . 1346088E409 

4.3053899E413 

Imaginary  Component 

6th  Order  0.000-0.100 

-2.E611943E+11 

3.1131480E+12 

-1.8348861E+13 

-2.4949310E+03  -2.4719206E+01 
1.9242115E+10  -6.8442981E+10 

2.4398E14E+06 

1.EE47730E+07 

-1.6171E97E+09 

G522(R) 


Real  Component 

6th  Order  0.001-0.100 

2.469086BE-01  8. 1E6471SE-02  -3.1310965E+02  -1 .7501311E+03  2. 1602490E+0B 

-2.66057ltE+06  9.6760956E+06 


Imaginary  Component 
8th  Order  0.001-0.010 

-5.0296804E+02  6.9944277E+05  -4.7773874E408  1.9210580E+11  -4.7861591E+13 

7.4464262E+16  -7.0205947E+17  3.6635682E+19  -8. 1050716E+20 

Ratio  0.010-0.100 

-3.0S06461E-04  3.7573808E-02  -1.6014618E+00  3. 1525211E-07  7.9630909E+02 

-5.0947594E-04  -7.9604799E+05  2.3725931E+07  -3.2364262E+08  2. 1545771E+09 

-5.6293092E+09 

Real  Component 

6th  Order  0.001-0.100 

9.6785892E-01  3.2169259E-01  -1 .2286658E+03  -6.8608196E+03  8.4780207E+05 

-1 . 0443582E+07  3 . 7986S23E*07 


Imaginary  Component 
8th  Order  0.001-0.010 

-7.1240104E+02  9.6940154E+05  -6.5262043E+08  2.6293123E+11  -6.5657068E+13 

1.0231241E+16  -9.6588315E+17  6.0443342E+19  -1 . 1166002E+21 
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Ratio  0.010-0.100 

-7.3279964E-04  8.4236082E-02 

-6.4107447E-04  -2.6362783E+06 
-6.0799782E+09 

G523(R) 

Real  Component 

6th  Order  0.000-0.100 

4.8908078E-01  9.7940015E-02 

-5.2961144E+06  1.9249510E+07 

Imaginary  Component 

8th  Order  0.000-0.010 

-1.0079854E+02  -3 . 9690570E+03 
■  -1.1930E63E+16  1 . 1354126E+17 

Ratio  0.010-0.100 

-4.6056133E-04  6.4394079E-02 

-4 . 4687389E-04  -6 . 1246793E+05 
-5.3094892E+09 

Real  Component 

6th  Order  0.000-0.100 

4.8908078E-01  9.794001BE-02 

-6.2961144E+06  1 .9249510E+07 

Imaginary  Component 

8th  Order  0.000-0.010 

-1.0079854E+02  -3.9690B70E+03 
-1 . 1930B63E+1B  1 . 13B4126E+17 

Ratio  0.010-0.100 

-4.60B6133E-04  6 . 4394079E-02 

-4 . 4687389E-04  -6 . 1246793E40B 
-S.3094692E409 


-3.2371602E+00  3.43B7684E-07  9. 13304B2E+02 

9.1863973E406  -1 .8340614E408  1.6197B94E-t-09 


-6.1B06941E402  -3.64B6986E-t-03  4.310173BE40B 


4.6810BB0E+07  -2.6387382E+10  7.36B6913E+12 

-B.8887B62E+18  1.2840182E+20 

-2.1168B82E+00  2.8B77B67E-07  8.4397927E+02 

1.8B33381E+07  -2.707129BE+08  1 . 9286742E+09 


-6.1B06941E+02  -3.6486986E+03  4.310173BE+08 


4.6810BB0E+07  -2.6387382E+10  7.36B6913E+12 

-B.8887B62E+18  1.2840182E+20 

-2.116S882E+00  2.8877B67E-07  8.4397927E+02 

1.8B33381E+07  -2.707129BE+08  1 . 9286742E+09 


Table  8.  Green’s  Function  Approximations  Summary 


Real  Component 

Imaginary  Component 

Function 

Approximation 

Range (m) 

Approximation 

Range (m) 

8th  Order 

Ratio  of  Polynomials 

0.001-0.010 

0.010-0.100 

6th  Order 

0.001-0.100 

GtM 

6th  Order 

0.001-0.100 

8th  Order 

Ratio  of  Polynomials 

0.001-0.010 

0.010-0.100 

G532(R) 

6th  Order 

0.000-0.100 

8th  Order 

Ratio  of  Polynomials 

0.000-0.010 

0.010-0.100 

8th  Order 

Ratio  of  Polynomials 

0.001-0.010 

0.010-0.100 

6th  Order 

0.001-0.100 

G^^(R) 

8th  Order 

Ratio  of  Polynomials 

0.000-0.010 

0.010-0.100 

6th  Order 

0.000-0.100 

G\zziR) 

6th  Order 

0.001-0.100 

8th  Order 

Ratio  of  Polynomials 

0.001-0.010 

0.010-0.100 

G523(R) 

6th  Order 

0.000-0.100 

8th  Order 

Ratio  of  Polynomials 

0.000-0.010 

0.010-0.100 

I[2i(R) 

8th  Order 

Ratio  of  Polynomials 

0.000-0.010 

0.010-0.100 

6th  Order 

0.000-0.100 

I*b2l(R) 

8th  Order 

Ratio  of  Polynomials 

0.000-0.010 

0.010-0.100 

6th  Order 

0.000-0.100 

8th  Order 

Ratio  of  Polynomials 

0.000-0.010 

0.010-0.100 

6th  Order 

0.000-0.100 

UR) 

8th  Order 

Ratio  of  Polynomials 

0.000-0.010 

0.010-0.100 

6th  Order 

0.000-0.100 

I*b3l(R) 

8th  Order 

Ratio  of  Polynomials 

0.000-0.010 

0.010-0.100 

6th  Order 

0.000-0.100 

8th  Order 

Ratio  of  Polynomials 

0.000-0.010 

0.010-0.100 

6th  Oruer 

0.000-0.100 
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Appendix  D.  Computer  Programs 


A  description  of  the  programs  and  subroutines  written,  modified,  and  used  during  this  research 
is  presented  in  this  appendix.  The  programs  written  by  Nazar  and  modified  by  the  author  were 
run  an  ELXSI  System  6400  computer  operating  under  the  ENBOS  13  operating  system.  The 
programs  written  by  the  author  were  run  on  a  20  Mhz  PC/AT  Clone,  Using  a  Microsoft  Version 
4.0  FORTRAN  compiler.  Copies  of  the  actual  code  may  be  obtained  by  contacting 

Major  Harry  Barksdale 
AFIT/ENG 

Wright-Patterson  AFB,  OH 
45433 

D.l  Programs  io  Numerically  Evaluate  the  Green’s  Functions 

GREENM  is  the  main  program  used  to  numerically  solve  the  integrals  of  all  the  Green’s 
functions  listed  in  Appendix  C.  All  the  material  parameters  of  the  antenna  are  set  in  GREENM. 
The  original  code  was  written  by  Nazar  [3]  and  only  evaluated  representative  examples  of  Green’s 
functions.  The  author  added  the  code  necessary  to  evaluate  remaining  Green’s  functions  examined 
in  this  thesis.  SUBROUT  contains  all  the  custom  written  subroutines  called  by  GREENM.  It  was 
also  written  by  Nazar  [3]  and  used  with  only  minor  modification.  EXTFUN  is  a  file  containing 
all  the  FORTRAN  implementation  of  the  integrands  for  the  various  Green’s  functions  and  all  non- 
IMSL  functions  used  in  GREENM.  It  also  was  originally  written  by  Nazar  and  modified  by  the 
author. 

D.2  Programs  Used  to  Approximate  the  Green’s  Functions 

LTSQR3  and  LTSQR8  are  custom  written  code  which  perform  a  least  squares  curve  fitting 
using  a  sixth  or  and  eighth  order  polynomial.  LTSQR82  performs  a  least  squares  curve  fit  using 
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an  eighth  order  polynomial  plus  two  terms  with  negative  powers.  The  result  is  a  curve  fit  using  a 
rational  function  composed  of  a  tenth  order  polynomial  divided  by  a  quadratic  polynomial.  All  three 
programs  call  two  subroutines  contained  in  SIMEQ  which  solve  a  system  of  simultaneous  linear 
equations.  The  routines  were  extracted  from  the  book  Numerical  Recipes  [21].  LTSQR6,  LTSQR8, 
and  LTSQR82  were  written  to  only  operate  on  the  data  from  one  component  of  a  Green’s  function  at 
a  time.  Lack  of  time  prevented  consolidation  of  all  the  routines.  Several  file  mamipulation  routines 
were  needed  to  break  out  the  input  file  components  and  combine  the  results  of  the  approximations. 
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